Introduction to the theory of left-handed media 
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This paper is devoted to investigating the physically interesting optical and electromagnetic 
properties, phenomena and effects of wave propagation in the negative refractive index materials, 
' which is often referred to as the left-handed media in the literature. This paper covers a wide 

range of subjects and related topics of left-handed media such as many mathematical treatment 
Mh . of fundamental effects (e.g., the reflection and the refraction laws on the interface between LH 

£ — ' and RH media, the group velocity and energy density in dispersive materials, the negative optical 

refractive index resulting from a moving regular medium, the reversal of Doppler effect in left- 
handed media, the reversal of Cerenkov radiation in left-handed media, the optical refractive index of 
massive particles and physical meanings of left-handed media, the anti-shielding effect and negative 
temperature in left-handed media, etc.), and their some applications to certain areas (e.g., three 
kinds of compact thin subwavelength cavity resonators (rectangular, cylindrical, spherical) made of 
left-handed media, the photon geometric phases due to helicity inversions inside a periodical fiber 
made of left-handed media, etc.). 

■ 



I. INTRODUCTION 



During the past 50 years, the design and fabrication of artificial materials attracted extensive attention in various 
scientific and technological areas [1]. For example, in the 1950s and 1960s artificial dielectrics were exploited for 
light-weight microwave antenna lenses, and in the 1980s and 1990s artificial chiral materials were investigated for 
microwave radar absorber applications [2]. During the last decades, a kind of material termed photonic crystals, 
which is patterned with a periodicity in dielectric constant and can therefore create a range of forbidden frequencies 
| called a photonic band gap focus considerable attention of many researchers [3] . Such dielectric structure of crystals 
04 1 offers the possibility of molding the flow of light inside media. Recently a number of interest was captured by the 
artificial electric and magnetic molecules [4] and artificial electric and magnetic materials, which, like many of the chiral 
materials, can exhibit positive or negative permittivity or permeability properties. More recently, a kind of artificial 
composite metamaterials (the so-called left-handed media) having a frequency band where the effective permittivity 
and the effective permeability are simultaneously negative receives particular and intensive attention of many authors 
G \ both experimentally and theoretically [5-14,49]. 

I ■ In the above, we present the briefly history of artificial materials, which was stated based on the point of view 
of pure technology. It is shown in what background (or environment) such a material with simultaneous negative 
permittivity and permeability arises. However, in what follows we will discuss this problem from the purely theoretical 
point of view. 

In 1930s, Dirac predicted a new particle (positive electron) based on his relativistic wave equation. In Dirac's 
theory, the energy eigenvalue of Fermions satisfies the equation E 2 = p 2 c 2 + m§c 4 with p and too denoting the 
momentum and rest mass of the Dirac particle. It follows that the energy eigenvalue E = ±y/p 2 c? + TOqC 4 . According 



to the convention of classical physics, the negative energy solution may be thought of as a one that has no physical 
meanings and therefore can be abandoned without any hesitation. But in quantum mechanics, the completeness of 
solutions is a very essential subject in solving wave equation. If we desert the negative energy solution, we cannot 
guarantee the completeness property of solutions of wave equation. So, the negative energy solution should also be 
of physical meanings and deserves consideration in quantum mechanics. Thus, the positive electron was predicted 
by Dirac. Likewise, in Maxwell electrodynamics, we may also experience such things. It is well known that in the 
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second-order Maxwellian equation, the relation between the squared of optical refractive index n and the permittivity 
e and permeability /i is n 2 = e/i. If both e and \x are positive, it follows that one can arrive at the following case of n: 

n = +y / eji, n = -y / qi, n = + y / {-e){-p), n = - A /(-e)(-/x), (1.1) 

all of which agree with the second-order Maxwellian equation. But for a time-harmonic electromagnetic wave, not all 
of them satisfies the first-order Maxwellian equation. In order to clarify this point, let us first consider the following 
two equations 

k x E = /i/iocjH, k x H = -ee ujE, (1.2) 

for a time-harmonic wave, where ^ an d eo stand for the electric permittivity and magnetic permeability in a free 
space, respectively. Here the wave vector k may be rewritten as k = n^k, where the unit vector k is so defined that 
k, E and H form a right-handed system. Thus according to Eq.(1.2), one can arrive at ^jj^k x E = H. It is readily 
verified that if fi > 0, then n is positive, and while if /i < 0, then n is negative. In the same fashion, from the equation 
k x H = — eenutE, one can obtain -2— k x H = — E. Thus the similar result can also be obtained: if n > 0, then e 
should be positive, and while if n < 0, then e should be negative. This, therefore, means that in Eq.(l.l), only the 
following two possibilities 

n = +v^u, n = -y/{-e)(-^) (1.3) 

may satisfy the first-order Maxwellian equation 1 . These two situations correspond to the two optical refractive indices 
of electromagnetic wave propagating inside right-handed media and left-handed media, respectively. Historically, in 
1967 2 , Veselago first considered this peculiar medium having negative simultaneously negative e and fi, and showed 
from Maxwellian equations that such media exhibit a negative index of refraction [16]. It follows from the Maxwell's 
curl equations that the phase velocity of light wave propagating inside this medium is pointed opposite to the direction 
of energy flow, that is, the Poynting vector and wave vector of electromagnetic wave would be antiparallel, i.e., the 
vector k, the electric field E and the magnetic field H form a left-handed system 3 ; thus Veselago referred to such 
materials as "left-handed" media, and correspondingly, the ordinary /regular medium in which k, E and H form a 
right-handed system may be termed the "right-handed" one. Other authors call this class of materials "negative- 
index media (NIM)" [17], "backward media (BWM)" [11], "double negative media (DNM)" [8] and Veselago's media. 
There exist a number of peculiar electromagnetic and optical properties, for instance, many dramatically different 
propagation characteristics stem from the sign change of the optical refractive index and phase velocity, including 
reversal of both the Doppler shift and Cerenkov radiation, anomalous refraction, amplification of evanescent waves 
[12], unusual photon tunneling [13], modified spontaneous emission rates and even reversals of radiation pressure to 
radiation tension [6] . 

In experiments, this artificial negative electric permittivity media may be obtained by using the array of long 
metallic wires (ALMWs) [18], which simulates the plasma behavior at microwave frequencies, and the artificial negative 
magnetic permeability media may be built up by using small resonant metallic particles, e.g., the split ring resonators 
(SRRs), with very high magnetic polarizability [19-21]. A combination of the two structures yields a left-handed 
medium. Recently, Shelby et al. reported their first experimental realization of this artificial composite medium, the 
permittivity and permeability of which have negative real parts [7]. One of the potential applications of negative 
refractive index materials is to fabricate the so-called "superlenses" (perfect lenses): specifically, a slab of such 
materials may has the power to focus all Fourier components of a 2D image, even those that do not propagate in a 
radiative manner [12,22]. 

In this paper, I will investigate theoretically the novel phenomena, effects and properties of wave propagation in 
left-handed media. 



1 Someone shows how the negative index of refraction arises: n = y / e/I — y (~ e )( — t 1 ) = \f~~^\/~P L = * 2 \/£M- We think that 
such a viewpoint may not be suitable for considering the existence of negative index of refraction. 

2 Note that, in the literature, some authors mentioned the wrong year when Veselago suggested the left-handed media. They 
claimed that Veselago proposed or introduced the concept of left-handed media in 1968 or 1964. On the contrary, the true 
history is as follows: Veselago's excellent paper was first published in Russian in July, 1967 [Usp. Fiz. Nauk 92, 517-526 
(1967)]. This original paper was translated into English by W.H. Furry and published again in 1968 in the journal of Sov. 
Phys. Usp. [16]. Unfortunately, Furry stated erroneously in his English translation that the original version of Veselago' work 
was first published in 1964. 

3 Note that by using the mirror reflection operation, a right-handed system may be changed into a left-handed one. It is thus 
clearly seen that the permittivity and permeability of the free vacuum in a mirror world may be negative numbers. 
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II. THE REFLECTION AND THE REFRACTION LAWS ON THE INTERFACE BETWEEN LH AND 

RH MEDIA 



The influence of the interface between LH and RH media on the radiation propagation is governed by the reflection 
and the refraction laws, which differ substantially from the laws in the conventional (right-handed) media. When a 
wave in a RH material hits an interface with a LH medium, it will have negative angle of refraction, i.e., the refracted 
wave will be on the same side of the normal as the incident wave, which are peculiar to the usual media. It is due to 
the fact that the refraction index is negative for LH media. Moreover, the refracted wave can be absent at all if the 
absolute values of the refraction indices of LH and RH media are equal [16,6] 



III. NEGATIVE PERMITTIVITY AND NEGATIVE PERMEABILITY 
A. Negative permittivity resulting from ALMWs structure 

A 2-layer periodic array of conducting straight wires will give rise to a negative permittivity. As observed, there 
is a single gap in the propagation up to a cutoff frequency for modes with the electric field polarized along the axis 
of the wire. Metals have a characteristic response to electromagnetic radiation due to the plasma resonance of the 
electron gas [18]. 

Here we consider the effective rest mass of electromagnetic wave in superconducting media (or electron plasma) 
and dispersive/absorptive materials. It is well known that the Lagrangian density of electromagnetics is 

< = 4^ F -^/¥ + wJ'A. (3.1) 
4 e c 

The equation of motion of an electron in superconducting media ( or electron plasma) acted upon by the external 
electromagnetic wave is m v = eE with E = — J^-A, where dot on the velocity v denotes the derivative of v with 
respect to time t and A stands for the magnetic vector potentials of the applied electromagnetic wave. One can 
therefore arrive at 

-| (m v + eA) = 0, (3.2) 

namely, the canonical momentum m c v + eA of the electron is conserved. Set m c v + eA = and then m v = — eA. 
So, the electric current density is of the form 

Ne 2 

J = New = -A (3.3) 

with TV being the electron number in per unit volume. It follows from (3.1) that the interaction term /j,q3 ■ A in 
Lagrangian density is rewritten as follows 

_ . UoNe 2 A 2 Ne 2 . 2 ,„ ^ 

Ai J-A = -^ A 2 = -A 2 , 3.4 

m c e m c c 2 

where use is made of the relation c 2 = j-^-- It is apparently seen that since electron current acts as the self- induced 

h 2 



eoAto 

charge current, the interaction term /ioJ ■ A is therefore transformed into the mass term —^0—A 2 of electromagnetic 
fields. Hence, we have 



=S£ = JtaL (3.5) 
h z e m c c 2 



and the effective mass squared 

2 H 2 ( Ne 



The relation between the frequency ui and wave vector k of lightwave propagating inside the dispersive medium may 
be jt = ~z with n being the optical refractive index. For highly absorptive media, the corresponding modification of 
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this method can also apply. So, for convenience, here we are not ready to deal with the case of absorptive materials. 

2 2 

This dispersion relation — yields 

up 1 

with the effective rest mass, m c g, of the light being so defined that m e ff, n and w together satisfy the following relation 

^-(1-nV. (3-8) 



?i 2 

This, therefore, implies that the particle velocity of light ( photons) is v = nc and the phase velocity of light is 

2 

v p = ^ = Since this relation is familiar to us, we thus show that the above calculation is self-consistent, at least 
for the de Broglie wave (particle). 

By using n 2 = e, it follows from Eq.(3.6) and (3.8) that 

uj 2 Ne 2 

e(w ) = l-% ul = ^- (3.9) 

with w p , N, e, m c and e being the electron plasma frequency, number of electrons per unit volume, electron charge, 
electron mass and electric permittivity in a vacuum, respectively. 

Generally speaking, when taking into consideration the absorption of media, the electric permittivity of ALMWs 
structures may be rewritten as follows 

oj 2 

e = l t — —. — r. (3.10) 

0J{W + l"f) 

Appropriate choice for the frequency of wave will yield negative permittivity. 

An alternative approach to deriving the expression for the permittivity of ALMWs structures can be seen in the 
Appendix A to this paper. 



B. Negative permeability resulting from SRRs structure 

The split ring resonator (SRR) can be used for manufacturing negative magnetic permeability media and hence 
designing left-handed materials. SRR is often formed by two coupled conducting rings printed on a dielectric slab of 
thickness with, say, 0.216 mm [23]. The structure of SRRs under consideration may be referred to the figures in the 
references [5,23]. The physical mechanism to obtain negative permeability may be as follows: a time- varying magnetic 
field applied parallel to the axis of the rings induces currents that, depending on the resonant properties of the unit, 
produce a magnetic field that may either oppose or enhance the incident field [5]. The associated magnetic field 
pattern from the SRR is dipolar. By having splits in the rings, the SRR unit can be made resonant at wavelengths 
much larger than the diameter of rings (the SRR particle size is about one tenth of that wavelength [23]). The purpose 
of the second split ring, inside and whose split is oriented opposite to the first, is to generate a large capacitance in 
the small gap region between the rings, lowering the resonant frequency considerably and concentrating the electric 
field [5]. 

In what follows we will derive the expression for the magnetic permeability of SRRs structures. Assume that the 
planes of SRRs are parallel to the x-y plane, and an external magnetic field B z = B| xt exp[ikj£]z acting on the SRRs. 
So, an electromotive force £ = —iwKr 2 B c ^ Lt (with nr 2 being the plane area of a SRR particle) is induced along the 
rings which is responsible for creating a current flow which produces a total magnetic moment in the particle. The 
slot between the rings acts as a distributed capacitance, which stores the same amount of charge (but of opposite 
sign) at both sides of the slot [23] . According to the charge conservation law, one can obtain 

itoCiVi - V ), 

iuC{V - Vi), (3.11) 



r &4> 
d 
r d(f> 



h = 



'at* 



-dt q ° 
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where <f> and C denote the angle displacement of 2D polar coordinate and the per unit length capacitance between 
the rings, and h and I the currents flowing on both rings (inner and outer), q; and q are the per unit length charge 
at the inner and outer rings, respectively. Note that on the right-handed sides of Eqs. (3.11), ^ can be replaced with 
iu>. It is assumed that in the ring the current can be considered to be homogeneous. This, therefore, implies that 
-p^Ii = It is apparent that (VI — V Q ) = 2iu (LI + i:r 2 B cxt ^, where L is the total inductance of the SRR. Thus, 



we have 



and consequently 



Thus the moment per unit volume is 



/ = 4TTiv 2 r C (LI + nr 2 B cxt ) (3.12) 



^r%C _ BCXt (313) 



1 - Airuj 2 LCr ' 



M = Nlnrl = ^( 4 ^:o-^o) jB ext > (314) 

1 — — 1 

where uj 2 — 4v l CrQ and N is the SRR unit number per unit volume. It can be easily verified 4 that the per unit length 
inductance of ring with radius r is ^rp. So, the total inductance of the SRR is Mo^o- It follows from Eq.(3.14) that 

Nnr 2 ^ 1 ^ 

M = _^B 0Xt . (3.15) 

1 — — T 

Let B lnt , H mt and H ext stand for the magnetic induction and field strength of internal and external fields, respectively 5 . 
Since all the currents (such as free current, polarization current and magnetization current) will contribute to magnetic 
induction B, while the magnetic filed H is produced only by the free current, we may have B lnt ^ _B cxt , H lnt — iJ cxt . 
Because of B int = MoMr# int = MoMr# cxt = ^B cxt , we can obtain B cxt = So, it follows that 



M= 2 "° . (3.16) 

1-4 Mr 



Insertion of M = ±!±=±B int yields 



In general, when taking account of the dissipation factor of SRR, the expression (3.17) may be rewritten as 

Flo 2 

Mr = 1 — — 2 , . r - (3.18) 

Appropriate choice for the frequency of wave will yield negative permeability. 

Readers may be referred to Appendix B for the derivation of the expression for the permeability of SRRs structures. 



4 In accordance with Ampere's circuital law, the magnetic induction B inside a long solenoid with N circular loops carrying a 
current I in the region remote from the ends is axial, uniform and equal to /Uo times the number of Ampere-turns per meter nl, 
i.e., B — /ion/, where n — -y- with I being the solenoid length. So, the magnetic flux is fionlTrrgl (here irr\ is the cross-section 
area of this long solenoid), and the electromotive force £ — —fxomvrol^. By the aid of £ = —L^, one can arrive at the total 
inductance of the solenoid, i.e., L = fj, nivrol. Thus the per unit length inductance of the solenoid is 2n r ni = 

5 B' nt refers to the field influenced by the moment. 
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In addition, someone may argue that the magnetization M should be M = J_ 1 jg cxt (rather than M = 
— ^ii— -B 1Dt ), where the relative permeability fi r is defined. Thus it follows from Eq.(3.15) that 



Nirr?" 1 



o^o z% 1 jU r 



l - ^0 Mr 



and in consequence 



which can be rewritten as 



which is also of the form analogous to (3.17). 



(3.19) 



1 ul-u? - (I + Nirr 2 Q )uj 2 + Nirrfa 2 

Mr i ^n 2 ^ 2 ^-(l + iVTrrg)^ 2 w 2 -(l + A^r 2 V 2 ' 1 ' 



C. Other design to obtain negative optical refractive index 

Pendry et al. used the array of long metallic wires (ALMWs) [18] that simulates the plasma behavior at microwave 
frequencies, and the split ring resonators (SRRs) that gives rise to the artificial negative magnetic permeability to 
realize the left-handed materials. But there are alternative approaches to the realization of such media. For example, 
Podolskiy et al. investigated the electromagnetic field disrtribution for the thin metal nanowires by using the discrete 
dipolc approximation, and considered the plasmon polariton modes in wires by using numerical simulation [24]. In 
the paper [24] , they described a material comprising pairs of nanowires parallel to each other and showed that such 
materials can have a negative refraction index in the near IR and the visible spectral range. Consider a thin layer of 
material, composed from pairs of nanowires parallel to each other. The length of individual nanowires is 2b±, their 
diameter is 262, and the distance between the nanowires in the pair is d. The needles are assumed to be embedded 
into host with dielectric constant. Podolskiy et al. considered the case of closely placed long nanowires so that 
b 2 <C d <C &i [24] . The incident wave propagates normal to the composite surface so that the electric field is parallel 
to the nanowires, while the magnetic field is perpendicular to the nanowire pairs. Podolskiy et al. found the effective 
dielectric constant and magnetic permeability for the 2D nanoneedle composite film are 

Ap d E Ap m H 

bib 2 d E 0i&2« H 

where p and mn are the surface metal concentration and the moment of the individual two-needle system, respectively. 
Note that the resonance position (and therefore the spectral range where the material refractive index is negative) is 
determined by parameters b\, 62 and d. By varying the these parameters the negative- refraction spectral range can 
be moved to the visible part of the spectrum [24]. 

It should be noted that the structure of a combination of two lattices (i.e., a lattice of infinitely long parallel wires 
and a lattice of the split ring resonators) is not isotropic. Such structures can be described with negative permittivity 
and permeability only if the propagation direction is orthogonal to the axes of wires [25] . More recently, Simovski and 
He presented a analytical model for a rectangular lattice of isotropic scatters with electric and magnetic resonances. 
Here each isotropic scatterer is formed by putting appropriately 6 £!-shaped perfectly conducting particles on the 
faces of a cubic unit cell [25] . They derived a self-consistent dispersion equation and used it to calculate correctly the 
effective permittivity and permeability in the frequency band where the lattice can be homogenized. The frequency 
range in which both the effective permittivity and permeability are negative corresponds to the mini-band of backward 
waves within the resonant band of the individual isotropic scatterer. For the detailed analysis of the electromagnetic 
properties of materials formed by the rectangular lattice of isotropic cubic unit cells of £1 particles, authors may be 
referred to Simovski and He's recently published paper [25] . 
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IV. WAVE PROPAGATION IN DISPERSIVE MATERIALS 



A. Group velocity and energy density in dispersive materials 



In a dispersive (isotropic) medium, the group velocity of an electromagnetic wave is 

duj 1 / 1 \ 



dfc 



where k = ^/ejl^€ a ^uj. In what follows we will calculate the energy flow velocity, 
Poynting vector may be given by 



|Ex H| = J—E 2 . 
V Wo 

Assume that the energy flow velocity equals the group velocity and we will obtain 

1 ( 1 \ 



_eeo_ £2 



(4.1) 

The magnitude of 
(4.2) 



w 



The right-handed side of Eq.(4.3) can be rewritten as 



(4.3) 



ee g2 



duj 



y 



The comparison between (4.3) and (4.4) shows that the energy density of electromagnetic system takes the form 



(4.4) 



(4.5) 



which can also be rewritten as 



w 



d d e 

-M + -(^)- 



(4.6) 



By using the expression H = . j^-E, -$-{u>ii)^-€qE' 2 is therefore 4-(u)i/)[IqH 2 . So, it follows from Eq.(4.6) that 



w = 



(4.7) 



which is just the electromagnetic energy density of dispersive linear materials. 
Note that according to the Kramers-Kronig relation 



Ree(u 



i + Ip 

7T 



+ °°^ld.', Ime M 



OO 



LU ' — LU 



-Ip 

7T 



[Ree(w') - 1] 



ui' — u> 



do/, 



(4.8) 



the negative refractive index materials is sure to be the dispersive (and also absorptive) media. So we should make use 
of Eq.(4.7) when calculating the energy density of electromagnetic materials. It can be easily seen that by substituting 
the expressions for the permittivity and permeability (3.10) and (3.18) of left-handed media into Eq.(4.7), the energy 
density is always positive. 

It is well known that in the "non-dispersive" medium, the electromagnetic energy density is w — i (E • D + H ■ B). 
How can we understand the relation between w and the electromagnetic potential energy (and kinetic energy) of 
dipoles? Does D contain the potential energy — P • E, or does it contain the kinetic energy of dipoles? 

Consider the following equation 



p + tu 2 p = eoUj^E 



(4.9) 
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that governs the electric dipole p in the presence of an electric field E = Eoexp [iwi] . The solution to Eq.(4.9) is 
p — po exp [ju>i\ , where 

— cj^ + sou;* 

Here = ^ , Wp = . Thus the kinetic and potential energies (purely mechanical energy) of the dipole oscillator 
(experience mechanical vibration) are written 

Ek = \^mx 2 = — ^muj 2 x 2 , 

E p = ^Kx 2 = ^muj 2 x 2 , (4.11) 
where x = -. So, the total mechanical energy of the dipole oscillator is 



£mccha tot = E k + E p = ~ (UJ 2 - LU 2 ) p 2 = (oj 2 - LU 2 ) p 2 . (4.12) 



The electric potential energy of the dipole oscillator is 



V = - P E = ^p 2 . (4.13) 

This therefore means that the relation between the mechanical energy and the electric potential energy of the dipole 
oscillator is 

V = -2£ mcchatot . (4.14) 
Thus the total energy of the dipole oscillator reads 

E to t= E mcchatot + V = -^ P E. (4.15) 

What is the relationship between E tot and ^pE in \ED with D = eoE+p! Now I will briefly discuss this problem: in 
fact, \pE in \ED is just the total mechanical energy, -E moc ha tot, of the dipole oscillator, while the electric potential 
energy has always been involved 6 in ^€qE 2 . This can also be understood as follows: when the incident wave with 
electric energy density \toE 2 in a vacuum propagates inside a medium, the electric energy density will be changed 
into \ED. Where does the energy difference \pE come? The answer lies in Eq.(4.15), which shows that the total 
energy (mechanical energy and electric potential energy) of the dipole oscillator is just the energy difference — \pE. 

It should be emphasized that when we take account of the electromagnetic energy density in the dispersive material, 
we should not neglect the terms | e w^E 2 + [Iqlu^H 2 . We think that in Ruppin's work [26], such a dispersive 
term should be added to his obtained expression for the electromagnetic energy density. 

B. Wave propagation in dispersive materials 

Now let us consider the wave propagation in a dispersive material, where the electromagnetic field equation that 
governs the wave propagation is 

V x £ + — (fiof/H) = 0, VxH-- (e e£) = 0. (4.16) 



Electric potential energy is involved in the electric field energy density. This may be discussed as follows: consider a 
self-interacting electric system, by using the Gaussian theorem V <j> — — we may obtain the total potential energy V = 

\ J p0dx = -±eo / 0V 2 0dx = -±eo / [V • (0V0) - (V0) 2 ] dx = ±e / £ 2 dx. 
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Here the electric and magnetic field strengths arc written in terms of quasi-monochromatic field components as follows 
[27] 

£ = cxp [— i(uj t — Pqz)] J cxp[— ia(t — /3'z)]E(x, uj + a)da, 

H = cxp [—i(uj t — Po z )] J cxp[—ia(t — /3'z)]H(x, uj + a)da, (4-17) 

where uj = loq + a, (3' = g^| Wo with ujq being the central frequency. Assuming the wave vector of electromagnetic 
wave is parallel to z-direction, it follows that 

d 

ifioGz x E + i[3'ac z x E + V x E — i^ wo/i( w o)H - i^a— (w^)H = 0, 

ouj 

d 

i/3 e z x H + i/3'ac z x H + V x H + ze woe(^o)E + ie a—(uje)~E = 0, (4.18) 

OUJ 

and consequently (multiplying the above two equations by H and E, respectively) 

d 

i/3 Ii ■ {c z x E) + i/3'atl ■ (e 2 x E) + H • (V x E) - i/i w M(wo)H 2 - in a—(Ljfi)H.' 2 = 0, 

OOJ 

d 

i/3 E • (e z x H) + i/3'aE ■ (e z x H) + E • (V x H) + ie u> e(w )'E, 2 + ie a^-(we)E 2 = 0. (4.19) 

OUJ 

Subtracting the second equation from the first one in (4.19), one can arrive at 

2i [foil ■ (e z x E) - e uj a e(uj )E 2 ] + V • (E x H) 



+ 2ia 



/3'H • (e 2 x E) - i Mo A (tJA1 )H 2 - I £o A (we ) E 2 



= 0. (4.20) 



It is readily verified that by using the relations /3 = V^"^ anc ^ I H| = y^^|E|, one can prove that the following 
relation 

/3 H • (e 2 x E) - e Q cu e(uj Q )-E 2 = (4.21) 
holds. Thus integrating the equation (4.20) and using J V • (E x H) dx = 0, we will obtain [27] 

i /_r s z d Z 



? C [^o^(^)H 2 + i £0 £ME 2 ] dz' 
Hence we obtain the expression for group velocity of wave propagation in the dispersive materials, i.e., 

It follows from (4.22) that the electromagnetic energy density w in dispersive media takes the form 



(4.22) 



(4.23) 



V. NEGATIVE OPTICAL REFRACTIVE INDEX RESULTING FROM A MOVING REGULAR MEDIUM 

Apart from the above artificial fabrications, does there exist another alternative to the negative index of refraction? 
In this section, I will demonstrate that a moving regular medium with a velocity in a suitable region (associated with 
its rest refractive index n) 7 possesses a negative index of refraction, which may be a physically interesting effect and 
might deserve further discussion. 



7 The rest refractive index is just the one measured by the observer fixed at the regular/ordinary medium. In what follows, 
we will derive the relativistic transformation for the optical refractive index (tensor) of a moving medium. 
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Suppose we have a regular anisotropic medium with the rest refractive index tensor n, which is moving relative to 
an inertial frame K with speed v in the arbitrary directions. The rest refractive index tensor h can be written as 
h = diag [m, ri2, 713] with rii > 0, i — 1,2,3. Thus the wave vector of a propagating electromagnetic wave with the 
frequency uj reads k = — {n\,n2,n^) measured by the observer fixed at this medium. In this sense, one can define a 
3-D vector n = (m, n 2 , n 3 ), and the wave vector k may be rewritten as k = n^. Now we analyze the phase tot — k • x 
of the above time-harmonic electromagnetic wave under the following Lorentz transformation 

x' = 7 (x-vi), f = 7 ( t -^E), (5.1) 

where (x', t') and (x, t) respectively denote the spacetime coordinates of the initial frame K and the system of moving 

_ i 

medium, the spatial origins of which coincide when t = t' = 0. Here the relativistic factor 7 = ^1 — . Thus 

by using the transformation (5.1), the phase u/t — k • x observed inside the moving medium may be rewritten as the 
following form by using the spacetime coordinates of K 



ojt 



-k-x = 7-(l-^)*'-7-("-J)-x', (5.2) 



the term on the right-handed side of which is just the expression for the wave phase in the initial frame K. Hence, 
the frequency uj' and the wave vector k' of the observed wave we measure in the initial frame K are given 

"' = 7-(l-^), k'= 7 ^-J), (5.3) 

respectively. 

To gain some insight into the meanings of the expression (5.3), let us consider the special case of a boost (of 
the Lorentz transformation (5.1)) in the ii-direction, in which the medium velocity relative to K along the positive 
cci-direction is v. In the meanwhile, we assume that the wave vector of the electromagnetic wave is also parallel to 
the positive aq-direction. If the rest refractive index of the medium in the aq-direction is n, its (moving) refractive 
index in the same direction measured by the observer fixed at the initial frame K is of the form 

n' = ^ = ^l (54 ) 

c 

which is a relativistic formula for the addition of "refractive indices" (— - provides an effective index of refraction). 
Note that here n' observed in the frame K may be negative. If, for example, when n > 1, the medium moves at 

c> v > - (5.5) 
n 

(and in consequence, k! > 0, uj' < 0); or when < n < 1, the medium velocity with respect to K satisfies 

c > v > nc (5.6) 

(and consequently, k' < 0, w' > 0), then n' is negative. The former case where k' > 0, to' < is of physical interest. 
In the paper [28] , it was shown that the photon in the negative refractive index medium behaves like its anti-particle, 
which, therefore, implies that we can describe the wave propagation in both left- and right- handed media in a unified 
way by using a complex vector field [28] . 

In the above, even though we have shown that in some certain velocity regions, the moving medium possesses a 
negative n' , such a moving medium cannot be surely viewed as a left-handed medium. For this point, it is necessary 
to take into account the problem as to whether the vector k, the electric field E and the magnetic field H of the 
electromagnetic wave form a left-handed system or not 8 . For simplicity and without the loss of generality, we choose 
the electric and magnetic fields of the electromagnetic wave in the medium system as 

E= (0,S 2 ,0), B = (0,0,B 3 ). (5.7) 



This requirement is the standard definition of a left-handed medium. Apparently, it is seen that the two concepts left-handed 
medium and negative refractive index medium is not completely equivalent to each other. 
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In the regular medium (right-handed medium with n > 0), the wave vector (along the positive ii-direction with the 
modulus k) and such E and B (or H) form a right-handed system 9 . In the frame of reference K, according to the 
Lorentz transformation, one can arrive at 

E' 2 = 7 (E 2 - vB 3 ) B' 3 = 7 (B 3 - ^E 2 ) . (5.8) 

For convenience, here we will think of the medium as an isotropic one, i.e., the refractive index n', permittivity e' 
and the permeability // of the moving medium are all the scalars rather than the tensors. First we consider the case 
of Eq.(5.6) where k' < 0, ui' > 0. In order to let the moving medium with a negative index of refraction be a real 
left-handed one, the wave vector (k' < 0), E' 2 and H' 3 of the electromagnetic wave measured in K should form a 
left-handed system. Because of k! < and B' 3 = fj,'^oH 3 with /j, being the magnetic permeability in a vacuum, if fi' 
is negative 10 , then only the case of E' 2 < and B' 3 > (and vice versa) will lead to the left-handed system 11 formed 
by the wave vector (k' < 0), E 2 and H' 3 . So, it follows from Eq.(5.8) that when the medium velocity relative to K is 
in the range 

(if El > c 2 Bl) , (5.9) 



(if El < c 2 Bl) , (5.10) 

then the wave vector, electric field and magnetic field measured in K will truly form a left-handed system. 

In conclusion, if the medium velocity satisfies both Eq.(5.6) and Eq.(5.9) or (5.10), then the moving medium will 
be a left-handed one observed from the observer fixed in the initial frame K. 

As far as the case of Eq.(5.5) where kl > 0, u>' < is concerned, it is also possible for the moving material to 
become a left-handed medium. But the case of u' < has no the practical counterpart (at least in the real situation 
where the artificial composite materials is designed). So, we would not further discuss the case of u)' < 0. 

Although the subject of the present note seems to be somewhat trivial, it is helpful for understanding the properties 
of wave propagation (such as causality problem, the problem of energy propagating outwards and backwards from 
source [29], etc.) in the left-handed medium. 



or 



JB 



3 v E 2 

E 2 B 3 



E 2 B 3 



VI. THE REVERSAL OF DOPPLER EFFECT IN LEFT-HANDED MEDIA 

As stated in the Introduction, the reversals of many optical and electromagnetic effects arise due to the negative 
index of refraction. In what follows we will consider these unusual phenomena. 

Let us assume that the medium is fixed at the initial frame K and the light source is moving inside the medium at 
velocity — v with respect to the rest frame K. It follows from Eq.(5.3) (i.e., w' = -fu; (l — that if here we consider 

a one-dimensional case in the regular (right-handed) medium with the refractive index n = 1, then the frequency we 

measure in the rest frame K is uj' = juj (l — ^) . Since the rclativistic factor 7 = ^1 — 2 , we may have 

which is the normal Doppler effect. However, if the medium is a left-handed material with n — — 1, then it follows 
from ui' = "fu (l — n ^ L ) that the frequency measured in K is u/ = juj (l + |) , which can be rewritten as 



c — v 

It is apparently seen that it is an abnormal Doppler effect compared with Eq.(6.1). 



(6.2) 



9 Thus if E'2 > is chosen to be positive, then B3 is also positive. In what follows, we will adopt this case (i.e., E2 > and 
B3 > 0) without the loss of generality. 

10 In the left-handed medium, the electric permittivity and the magnetic permeability are simultaneously negative. 

11 Since k' < 0, in order to form a left-handed system, the signs of E' 2 and H' 3 should not be opposite to each other. So, the 
signs of E' 2 and B 3 should be opposite due to fif < 0. 
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VII. THE REVERSAL OF CERENKOV RADIATION IN LEFT-HANDED MEDIA 



Electromagnetic radiation, usually bluish light, emitted by a beam of high-energy charged particles passing through 
a transparent medium at s speed greater than the speed of light in that medium. This effect is similar to that of a 
sonic boom when an object moves faster than the speed of sound. In this case the radiation is a shock wave set up in 
the electromagnetic field. 

In the theory of Cerenkov radiation, the Fourier component of the magnetic fields produced by the moving charged 
particle is 

B u = - A 3 p ; sin 05 COS0), 7.1 

4:TT€ a C 6 R V V c ) 

where 0is the angle between the particle speed v and the direction of radiation. Let us first assume that the refractive 
index n > 0. It is well known that when the charged particle velocity v satisfies v > ^, then in the direction of 
cos0 = — , the magnetic field B u is nonvanishing. If the refractive index n < 0, then the radiation direction may 
change, i.e., 9 — > 9 + ir. Additionally, the reversal of the direction of will also arise, which results from the 
coefficient 4 ^"^ 3 of (because of the coefficient involving n). 



VIII. OPTICAL REFRACTIVE INDEX OF MASSIVE PARTICLES AND PHYSICAL MEANINGS OF 

LEFT-HANDED MEDIA 

A number of novel electromagnetic and optical properties in left-handed media result from the negative index 
of refraction. But what is the physical meanings of the negative index of refraction? Although in the literature 
many researchers have investigated the negative refractive index of left-handed media by a variety of means of applied 
electromagnetism, classical optics, materials science as well as condensed matter physics, less attention than it deserves 
is paid to the physical meanings of negative refractive index. In this section we will study this fundamental problem 
from the purely physical point of view. 

Before we consider the negative refractive index of left-handed media, we first deal with the "optical refractive 
index" problem of de Broglic wave, the results of which will be helpful in discussing the former problem. 

For the case of de Broglie particle moving in a force field, we can also take into account its "optical refractive index" 
n. According to the Einstein- de Broglie relation, the dispersion relation of the de Broglie particle with rest mass m 
in a scalar potential field V (x) agrees with 



c 2 c 2 + 



(8.1) 



where <j> is defined to be <fi — and k, h and c stand for the wave vector, Plank constant and speed of light in a 
vacuum, respectively. It follows from Eq.(8.1) that 



2- 

LU 



2„2 



k'c 



(8.2) 



which yields 



m^c 4 



2^ + 



(i) 



2 ■ 



(8.3) 



Compared Eq.(8.3) with the dispersion relation ^ 



one can arrive at 

2 



n 2 = l 



h 2 ^ 



.0 



(8.4) 



which is the square of "optical refractive index" of de Broglie wave in the presence of a potential field 4>. 

It is of physical interest to discuss the Fermat's principle of de Broglie wave in a potential field V. One can readily 
verify that the variation 5 J cn 2 dt = can serve as the mathematical expression for the Fermat's principle of de Broglie 
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particle (for the derivation of this expression, the readers may be referred to the the paper [28]). If, for example, by 
using the weak-field and low-motion approximation where 



< 1, < 1 and (fy in Eq.(8.4) can therefore be 



ignored, 1 — j-^jj — 72-, then n 2 is approximately equal to 



n 2 ~ f imu 2 - V). (8.5) 
me 2 V2 / V ; 

Thus in the non-relativistic case, Fermat's principle, S J cn 2 dt — 0, is equivalent to the following well-known action 
principle S J L(x, v)dt = (differing only by a constant coefficient ^5), where the function L(x, v) = imv 2 — V(x) 
denotes the Lagrangian of massive particle in a potential field V(x). Hence the above formulation for the "optical 
refractive index" of the de Broglie particle inside a potential field is said to be self-consistent. 

According to Eq.(8.4), there are two square roots (i.e., positive and negative roots, n + , n_) for the refractive index 
of de Broglie wave. It is well known that a subatomic particle has its corresponding antiparticle that has the same 
mass as it but opposite values of some other property or properties. Indeed, it is verified that the particle and its 
antiparticle possess these two square roots, respectively. If the refractive index of particle is positive, then that of its 
antiparticle will acquire a minus sign, which may be seen by utilizing the charge conjugation transformation: uu — > — u, 
<p — > —<j). If, for example, we take into consideration the refractive index of electron and positron by regarding the 
positron as the negative energy solution to Dirac's equation, then we can classify the solutions of Dirac's equation 
into the following four cases: 

(i) E = E + , h = hk, n = n + ; 

(ii) E = E + , h=—hk, n = n + ; 

(iii) E = E-, h = hk, n = n_; 

(iv) E = E-, h=—hk, n = n-, 

where h, E + and respectively represent the helicity, positive and negative energy corresponding to the electron 
and positron. It is apparent that the relationship between the cases (i) and (iv) is just the charge conjugation 
transformation. This transformation also relates the case (ii) to (iii). 

Compare the case of wave propagation in left-handed media with the results presented above, we can conclude that 
the negative optical refractive index in left-handed media corresponds to the antiparticles of photons (the so-called 
antiphotons) . However, as is well known, there exist no antiphotons in free space. The theoretical reason for this is that 
the four-dimensional vector potentials A n with jj, = 0, 1, 2, 3 are always taking the real numbers. But in a dispersive 
and absorptive medium, as an effective medium theory [30], the vector potentials of which may probably take the 
complex numbers. Such complex vector field theory has been considered previously [31]. Here the Lagrangian density 
is written C = — ^F* u F^ Vl where the electromagnetic field tensors — d^A v — d u A Ul F* v = d^A* v — <9„A* with A* 
being the complex conjugation of A^. The complex four-dimensional vector potentials A^ and A* characterize the 
propagating behavior of both photons and antiphotons. 

Now a problem left to us is that does the complex vector field theory apply well in investigating the light propagation 
inside the negative refractive index medium. Detailed analysis will show that this theory is truly applicable to the 
consideration of optical index of refraction and wave propagation in left-handed media. This will be proved in what 
follows by using another mathematical form 

~M = VxM, V-M = (8.6) 

for the Maxwellian equations, where M = + i^JIH, c n = e and [i are absolute dielectric constant and 

magnetic conductivity of the medium, respectively. It can be easily verified that the Maxwellian equations can be 
rewritten as Eq.(8.6). 

If e > and [i > are assumed, then the light propagation inside the isotropic linear right-handed medium (regular 
medium) is governed by Eq.(8.6). For the time-harmonic electromagnetic wave in the right-handed medium, Eq.(8.6) 
can be rewritten 

i-M = -kxM (8.7) 

with cu being the frequency of electromagnetic wave. 

Now let us take the complex conjugation of the two sides of Eq.(8.7), and the result is of the form 

i— M* = k x M*, (8.8) 
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where the field M* is expressed in terms of both electric and magnetic fields, E* and H*, i.e., M* = ^eE* — i^/jlH* . 
Thus, we have 

i — [ViE* - i^/MH*] = k x [>/eE* - *VmH*]. (8.9) 
Multiplying the two sides of Eq.(8.9) by the imaginary unit i, one can arrive at 

i — [V~eE* - iV^JlH*} = k x [V~eE* - iy/=jiH?], (8.10) 

where the magnitude of c„ (i.e., l/y/eji) docs not change but it can be rewritten as a new form c„ = e)(— fi). 

According to Eq.(8.10), one can obtain the following relations 

kxE* =w(-/i)H*, k x H* = — w(— e)E*. (8.11) 

Note that here e > 0, /i > 0. Thus it follows from Eq.(8.11) that the wave vector k, the electric field E* and the 
magnetic field H* of the electromagnetic wave in this medium, through which the fields E* and H* propagate, form a 
left-handed system. This, therefore, means that here the Poynting vector S is pointed opposite to the direction of the 
wave vector of electromagnetic wave. For this reason, it is expected that the reversal of wave vector, instantaneous 
helicity inversion [32] and anomalous refraction will take place when an incident lightwave from the right-handed 
medium travels to the interfaces between left- and right- handed media. Since the direction of phase velocity and 
energy flow of lightwave propagating in the left-handed medium would be antiparallel, the change of wave vector k will 
truly occur during the light propagation through the interfaces. The inversion of k means that its magnitude acquires 
a minus sign, namely, the optical index of refraction n becomes negative since the magnitude of k is defined to be 
k = nuj/c. In view of the above discussion, it is concluded that Eq.(8.10) governs the propagation of time-harmonic 
planar electromagnetic wave, where the electromagnetic fields are characterized by E* and H*. 

Thus by making use of the complex vector field theory we obtain a formulation which can treat the wave propagation 
in both right- and left- handed media. One of the advantages of the present formulation is that the minus sign of 
negative refractive index can be placed into the theoretical calculations by using the above unified approach rather 
than by hand. Such unified approach will be useful in dealing with the electromagnetic wave propagation and photonic 
band gap structures in the so-called left-handed medium photonic crystals, which are artificial materials patterned 
with a periodicity in negative and positive indices of refraction. 



IX. ANTI-SHIELDING EFFECT AND NEGATIVE TEMPERATURE IN LEFT-HANDED MEDIA 

In linear media the electric and magnetic polarizations (per unit volume) is respectively 

P = (e-l)e E, M=0u-1)H. (9.1) 

The potential energy density of electric and magnetic dipoles in electromagnetic fields reads 

[/ e = -P-E = -( e -l)e E 2 , U m = -;i M • H = —(fi - l)/i H 2 (9.2) 

with no being the permeability of free space. It is apparent that in left-handed media where both e and [i are negative, 
the electric and magnetic polarizations, P and M, are opposite to the electric and magnetic fields, E and H, thus in 
general, the potential energy density, both U e and U m , are positive, which differs from the cases in the conventional 
media (i.e., the right-handed media). This, therefore, means that left-handed media possess the anti-shiclding effect, 
which results in many anomalous electromagnetic and optical behaviors [6], as is mentioned above. 

May the electric and magnetic dipole systems be brought into the state of negative absolute temperature when 
the electromagnetic wave is propagating in left-handed media? We show that under some restricted conditions, these 
dipole systems may possess a negative temperature indeed. The detailed analysis is given in what follows: 

According to Ramsey's claim [33] that the essential requirement for a thermodynamical system to be capable of a 
negative temperature is that the elements of the thermodynamical system must be in thermodynamical equilibrium 
among themselves in order that the system can be described by a temperature at all. For the left-handed media, 
this may be understood in the sense: the electric and magnetic dipole moments produced by ALMWs and SRRs 
should be coupled respectively to each other, and this dipole-dipole interactions must be strong and therefore the 
thermodynamical equilibriums between themselves may be brought about in a short time, which is characterized 
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by the relaxation time t\. But here both ALMWs and SRRs are macroscopic elements of left-handed medium 
(e.g., the length scale of SRRs is mm), which implies that there exists a problem as to whether the concept of 
thermodynamical equilibrium is applicable to these macroscopic elements or not. This problem, however, is not the 
subject of this section, and it will be published in a separate paper; furthermore, if, for example, we can find a kind 
of molecular- type SRRs (i.e., the C- or il- like [34] curved organic molecule that possesses the delocalized n bond and 
can therefore provide a molecular electric current induced by magnetic field of lightwave) for manufacturing the left- 
handed media, then the dipole system of this media can be appropriately described by a temperature [33] . So in our 
tentative consideration in this paper, we can study the anti-shielding effect and negative temperature of left-handed 
media from the phenomenological point of view, ignoring the detailed information about the elemental atoms which 
constitute the left-handed media. 

In order to generate a state of negative absolute temperature, the following requirement for a thermodynamical 
system should also be satisfied: the system must be thermally isolated from all other systems, namely, the thermal 
equilibrium time (ti) among the elements of the negative-temperature system must be short compared to the time 
during which appreciable energy is lost to or gained from other systems [33] . In left-handed media, the dipole- lattice 
interaction causes heat to flow from the dipole system (negative temperature) to the lattice system (positive tempera- 
ture). If the relaxation time (T2) during which this interaction brings the dipole-lattice system into thermodynamical 
equilibrium is much greater than t± , then it may be said that the dipole system is thermally isolated from the lattice 
system, and the dipole system is therefore capable of a negative temperature. In addition to this condition, there 
exists another restriction, i.e., the state of negative temperature should be built up within half period of electromag- 
netic wave propagating inside the left-handed medium, i.e., the relaxation time n should be much less than j with / 
being the frequency of electromagnetic wave, which means that this thermodynamical process caused by dipole-wave 
interaction can be considered a quasi-stationary process and in consequence the state of negative temperature of 
dipole system due to dipole-dipole interaction may adiabatically achieve establishment in every instantaneous interval 
during one period of electromagnetic wave. If this condition is not satisfied, say, the relaxation time n is several 
times as long as j, then the negative temperature of dipole systems cannot be achieved since the thermodynamical 
equilibrium among dipolcs induced by the oscillating electromagnetic fields E and H may lose its possibility, namely, 
the concept of thermodynamical equilibrium is not applicable to this non-stationary process; in fact, it should be 
treated by using non-equilibrium statistical mechanics. In view of above discussions, we conclude that the principal 
conditions for the electric and magnetic dipole systems of left-handed media possessing a negative temperature may 
be ascribed to the following two essential requirements 

ti < t 2 , n < j. (9.3) 

Generally speaking, the first restriction, n <t 2 , is readily achieved, since in electromagnetic media the dipole-dipole 
interaction are often much stronger than the dipole-lattice interaction. For example, in some nuclear spin systems 
such as the pure crystal of LiF where the negative temperature was first realized in experiments, the relaxation time 
n of spin-spin process is approximately the period of the Larmor precession of one nucleus in the field of its neighbor 
and is of the order of 10~ 5 second while the relaxation time that depends strongly upon the interaction between the 
spin system and the crystal lattice is several minutes [33]. In the recent work of designing the structures of more 
effective left-handed media, use is made of the SRRs of high electromagnetic polarizability [5]. This leads to the 
very strong dipole-dipole interaction and may achieve the thermodynamical equilibrium of dipole system with shorter 
relaxation time. If this interaction is rather strong and therefore agrees with the second restriction, t\ <C j with / 
being in the region of GHz [18-20], then the dipole systems of left-handed media may be said to be in the state of 
negative absolute temperature. 

Since the most conspicuous feature of left-handed media may be the electromagnetic anti-shielding effect, the 
previous considerations concerning the superluminal light propagation in the instantly varying electric fields based 
on the viewpoint of effective rest mass can also be applied to the incident electromagnetic wave propagating inside 
left-handed media. As is stated, in left-handed media many electromagnetic effects such as Doppler shift, Cherenkov 
radiation and wave refraction are inverted compared with those in ordinary right-handed media. In the similar fashion, 
it is believed that this superluminal light propagation may also probably occur in left-handed media. Moreover, in 
some electromagnetic materials, there exists another superluminal pulse propagation that also results from the anti- 
shielding effect [35,36] (but rather than from the statistical fluctuations in negative temperature), since in this media 
the designed equivalent permittivity and permeability based on the so-called two-time-derivative Lorentz model are 
simultaneously smaller than the values in free space [35,36]. Ziolkowski has studied the superluminal pulse propagation 
and consequent superluminal information exchange in this media, and demonstrated that they do not violate the 
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principle of causality [1]. We suggest that all these potential supcrluminality effects of light propagation in anti- 
shielding media should be further investigated, particularly in experiments. 



X. FREQUENCY-INDEPENDENT EFFECTIVE REST MASS OF PHOTONS IN THE 2TDLM MODEL 

A. Extracting frequency-independent effective mass from refractive index squared 

In this framework of definition of frequency-independent effective rest mass, the photon should be regarded as acted 
upon by a hypothetical force field, namely, it behaves like a massive de Broglie particle moving in a force field. It has 
been verified previously, for the case of de Broglie particle in a force field, one can also consider its "optical refractive 
index" n. 

In what follows, an approach to frequency- independent effective mass extracted from the optical refractive index 
squared n 2 (w) is presented. It is assumed that n 2 (u>) can be rewritten as the following series expansions 



- 2 H = E^< ( 10 -D 

k 

then one can arrive at n 2 (uj)lo 2 — J2k ~£*=z> namely, 

n 2 (uj) lo 2 = ... + — - + a_ 2 + a_iUO + a ui 2 + a +1 w 3 + ... . (10-2) 

UJ 

Compared (10.2) with (8.4), it is apparent that the frequency-independent effective rest mass squared m 2 ^ of the 
photon can be read off from n 2 (oj), namely, the constant term on the right-handed side on (10.2) is related only to 
m 2 cS , i.e., 



«-2 = -^- (10.3) 

h 2 



Thus it follows from Eq.(8.4) that after we extract the term co "| t ' , the retained terms 1Z in n 2 (oj) — 1 belongs to 
—2^ + (jL^J . If the effective potential <fi can be rewritten in a form <f> = J2k §kU k with <pk being ^-independent, then 

one can obtain the expansion coefficients (f>k of <p. It should be noted that here the terms in —2^ + (j^j that takes 
the form co "| t - has already been involved in We will clarify further this point in what follows. It is apparently 
seen that by using the form of series expansion of <f>, the terms —2^ and (j^J in Eq.(8.4) can be rewritten 



-2* = -2*± -2 £ (10.4) 

and 



u> u> 



'" + 2Efc 2 >0 ^ fc )+E^ 2fc ' 2 + 2 E ^ k+l - 2 - (io.5) 

' k^O k^-Lk>l 

Note that the coefficient factors of uj~ 2 , i.e., —2<f>^i in Eq.(10.4) and 4"o + 2 Tl l k>o < t ) k4>-k in Eq.(10.5) contributes 
to the photon effective rest mass, that is, it is difficult to distinguish the effect of these coefficient factors from the 
effective rest mass. So we can attribute the effective rest mass term to the coefficient factors of w~ 2 , namely, a_ 2 
contains — 20_i + <j) 2 + 2 J2k>o ^k^-k- After the term co J| t ' is extracted from n 2 (uj) — 1, the retained term 1Z equals 
the summation of those terms in Eq.(10.4) and (10.5) that does not take the form c< ^| t- , i.e., 

K = -2j2^ k - 1 +J2^ 2k ~ 2 + 2 E (10-6) 

kyt-l k^tO k^-l,k>l 

To conclude, both the photon effective rest mass m c g and the effective potential <p are defined in the above for- 
mulation. In a free space, since n 2 — 1 = and the vacuum contributes no effective mass to photons. However, in a 
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dispersive medium, if n 2 (uj) — 1 is regarded as the contribution of the interaction between light and medium, then the 
term co "| t - in n 2 (uj) — 1 will result in an ^-independent photon effective rest mass. Other terms retained in n 2 (uj) — 1 
can be ascribed to the effective potential cj>, which is so defined that its expansion coefficients cf) k satisfy Eq.(10.6). 

Generally speaking, the optical refractive index n 2 (uj) is often of complicated form, particularly for the metama- 
terials. So, we cannot extract frequency-independent m 2 ff immediately from n 2 (to). In the following, this problem is 
resolved by calculating the limit value of n 2 (uj) uj 2 , where uj tends to oo or zero. We first consider the large frequency 
behavior of n 2 (uj)uj 2 . It follows from (10.2) that n 2 (uj)uj 2 at high frequencies (i.e., uj — > oo) is of the form 

n 2 (uj) u?^a-2 + fl-iw + a a uj 2 + a + iw 3 + a +2 w 4 + ... . (10-7) 

This, therefore, means that n 2 (uj) uj 2 at high frequencies (uj — > oo) can be split into the following two terms 

n 2 (uj) uj 2 ^> (constant term) + (divergent term) . (10.8) 

So, we can obtain m 2 s from the following formula 

2 4 

lim \n 2 (uj)uj 2 - (divergent term)] = - m ° sC , (10.9) 

Note that the constant term in n 2 (uj) uj 2 can also be extracted by considering the low-frequency (zero-frequency) 
behavior of n 2 (uj). It is certain that in this case the obtained constant term is the same as (10.9). Apparently 
n 2 (uj) uj 2 at low-frequencies (i.e., uj — ► 0) behaves like 

n 2 (uj)uj 2 ^... + ^ + ^- + ^+a_ 2 , (10.10) 

UJ 6 UJ Z UJ 

where the sum term (... + + + ^-) is divergent if uj approaches zero. Thus, subtraction of the divergent term 
from n 2 (uj) uj 2 yields 

lirn [n 2 (uj) uj 2 - (divergent term)] = - . (10.11) 

In the next section, we will calculate the frequency-independent effective rest mass of photons in the two time 
derivative Lorentz material (2TDLM) model [35,36] by making use of the method presented above. 



B. Frequency-independent effective rest mass of photons in 2TDLM model 

One of the metamaterials models proposed by Ziolkowski and Auzanneau is called the two time derivative Lorentz 
material (2TDLM) model [35,36], which is a generalization of the standard Lorentz material model and encompasses 
the permittivity and permeability material responses experimentally obtained, which may be seen, for example, in 
the reference of Smith et al. [5] . According to the definition of the 2TDLM model, the frequency-domain electric and 
magnetic susceptibilities are given [1] 

= (<) 2 X^^<x^ 2 x; m M = K m ) 2 x^ + ^"X""^XT , 12) 

-uj 2 + iuT* + (ulf ' -uj 2 + iujT™ + (uj™) 2 

where Xa' ro i X^'™ an d X^'™ represent, respectively, the coupling of the electric (magnetic) field and its first and second 
time derivatives to the local electric (magnetic) dipole motions. uj p , T e , T m and uj can be viewed as the plasma 
frequency, damping frequency and resonance frequency of the electric (magnetic) dipole oscillators, respectively. So, 
the refractive index squared in the 2TDLM model reads 

n 2 (w) = 1 + x e M + x m M + x e M x m M • (10.13) 

In order to obtain m 2 s from n 2 (uj), according to the formulation in the previous section, the high-frequency (uj — > oo) 
behavior of x e ( w ) ^ 2 should first be taken into consideration and the result is 

X e (uj)uj 2 ^ -iujuj; X e p + u 2 x e r (10.14) 
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Note that — iu)ui^x% + w X 7 is the divergent term of x e ( w ) w a t large frequencies. Subtracting the divergent term 
from x e (w) uj 2 , we consider again the high-frequency behavior (uj — > 00) and obtain another divergent term iw 3 r e x 7 , 
i.e., 

X e (w) - (-iww;x| + ^ 2 x 7 ) 
_ ^ 2 (^g) 2 Xq - ^ 2 ^r e x| + io^g (u;g) 2 - iu z Y e x% ~ u 2 K) 2 X 7 

Thus, after subtracting all the divergent terms from x e (w) w 2 , we obtain 

jim [x e m^ 2 - + ^ 2 x 7 ) - ^ 3 r e x;] 



(10.15) 



(c,;) 2 xl ^ e x% - H) 2 x 7 + (r e ) 2 x^ 



(10.16) 



In the same fashion, we obtain 



hm [x m (w) c 2 - (-iww?x? + " 2 x7) - ^ 3 r m x7] 
( Wp ") 2 x™ - <r m x™ - K m ) 2 x7 + ( r " 1 ) 2 x 



(10.17) 



for the magnetic susceptibility x" 1 (w). In what follows we continue to consider the term x e (w) x" 1 ( w ) on the right- 
handed side of (10.13). In the similar manner, the infinite-frequency limit of x e ( w ) X™ ( w ) ^ 2 is given as follows 



hm x e (w) x m M ^ = -x™ K) 2 x„ - ^Tx% Hf x 7 + (r e ) 2 x 7 



-x 



J pJ - ^pi X/3 - 1^0 j X 7 + l 1 j X 7 



(10.18) 



Hence, insertion of (10.16)-(10.18) into (10.9) yields 



'eff- 



= (1 + x7) Hfxl - ^ e x e P K) 2 x 7 + (r e ) 2 x\ 



i + x 7 



( W P J X Q - W P 1 X/3 - (Vo J X 7 + I 1 J X 7 



(10.19) 



Thus the w-independent effective rest mass squared m 2 ff is extracted from n 2 (u>). Apparently, m 2 ff in Eq.(10.19) 

2 4 

does not depend on u>. The remainder on the right-handed side of Eq. (8.4) is n 2 (to) — 1 + 2 ; which is equal to 
(w) — ^w - ^ or simplicity, we set 



'eff- 



7i /. 



(10.20) 



and then from the equation 

g) 2 -2^A H , (10.21) 

one can readily obtain the expression for the hypothetical potential field (f> (lu) is (j> (lu) = uj 1 — yjl + A (uj) , or 

(10.22) 



V (w) = TlUJ 1 - y/l + A(uj) 

where <j> = ^. 

Many metamaterials which have complicated expressions for electric and magnetic susceptibilities such as (10.12) 
are generally rarely seen in nature and often arises from the artificial manufactures. In the following, we will consider 
the photon effective rest mass in one of the artificial composite metamaterials, i.e., the left-handed medium. In 
general, the dielectric parameter e (uj) and magnetic permeability (lu) in the isotropic homogeneous left-handed 
medium are of the form 
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*(«,) = 1- 2 r (10-23) 

with the plasma frequency u> v and the magnetic resonance frequency luq being in the GHz region. 7 and T stand 
for the electric and magnetic damping parameters, respectively. Compared (10.23) with the electric and magnetic 
susceptibilities (10.12), one can arrive at 

Xa = 1. w p = w p> x% = x% = 0, = 0, r e = -7, 

x™ = -f, x™ = x% = o, ^ = ^, r m = -r. (10.24) 

So, it is readily verified with the help of (10.19) that the w-independent effective rest mass of photons inside left-handed 
media is written as follows 



m off f 

h 2 



(1-F)^ + F(^_r 2 ). (10.25) 



It follows that Eq.(10.25) is a restriction on the electromagnetic parameters F,uj p ,oj n ,T and F in the electric 
permittivity and magnetic permeability (10.23), since m 2 s > 0. Insertion of the experimentally chosen values of the 
electromagnetic parameters in the literature into Eq. (10.25) shows that this restriction condition is satisfied. For 
instance, in Ruppin's work [37], F, w p ,wo,7,r and F were chosen F — 0.56, ui p — 10.0GHz, loq = 4.0GHz, 7 = 0.03o> p , 
T = 0.03a>cb which indicates that the inequality m 2 s > is satisfied. 



C. Some applications of frequency-independent effective rest mass of photons 

In this section we consider the applications of the concept of frequency-independent effective rest mass of photons 
to some optical and electromagnetic phenomena and effects. 

(i) As was stated above, since the frequency-independent effective rest mass of photons in electromagnetic media 
is related close to the coupling parameters of the electric (magnetic) fields to the electric (magnetic) dipole systems, 
and the damping and resonance frequency of the electric (magnetic) dipole oscillators, it contains the information 
on the interaction between light and materials. It is shown that in certain frequency ranges this effective rest mass 
governs the propagation behavior of electromagnetic wave in media, and by using this concept one can therefore 
consider the wave propagation somewhat conveniently. For example, in the left-handed media whose permittivity and 



permeability are expressed by (10.23), the potential field V(uj) approaches —Ti^J F(uj 2 — u 2 ), which is constant, as the 
frequency tends to zero. This, therefore, means that at low frequencies the light propagation is governed mainly by 
the frequency-independent effective rest mass of photons. 

Likewise, in the case of the potential field V(co) being vanishing where the frequency squared of electromagnetic 
wave is 



Q 2 = , (10.26) 

the photon can be regarded as a "free" massive particle with the rest mass being expressed in (10.25). In the frequency 
region near fl, sine the photon is a "quasi-free" particle, the light propagation behavior and properties is rather simple 
compared with other frequency ranges in which V(w) is nonvanishing or large. For the left-handed media, the relevant 
frequency parameters luq and ui p of which is about GHZ, the frequency fl in Eq. (10.26) that leads to V(Q.) = is 
just in the frequency band (GHZ) where the permittivity and permeability are simultaneously negative. It follows 
that when the effects and phenomena associated with negative index of refraction occur in left-handed media, the 
wave propagation (e.g., scattering properties) and dispersion properties of the light are not as complicated as those 
in other frequency ranges, since in the negative-index frequency band, the optical refractive index squared simulates 
nearly the plasma behavior, particularly in the region near O. 

(ii) Historically, possibility of a light pulse with speed greater than that (c) in a vacuum has been extensively 
investigated by many authors [1,38-42,45]. Ziolkowski has studied the superluminal pulse propagation and consequent 
superluminal information exchange in the 2TDLM media, and demonstrated that they do not violate the principle of 
causality [1]. He showed that in the 2TDLM model, both the phase and group velocities, i.e., 

lim v(lu) ~ -. — - — r (10.27) 
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of light at high frequencies will exceed the speed (c) of light in a vacuum, so long as —2 < X7 < 0, where \i denotes 
the coupling coefficients of the second time derivatives of electric (magnetic) fields to the local electric (magnetic) 
dipole motions [1]. 

Although in the 2TDLM model the medium will exhibit a possibility of superluminal speeds of wave propagation, 
the photon velocity (rather than group velocity of light) is by no means larger than c. This may be explained as 
follows: according to Eq.(10.22), the kinetic energy of a photon is £ = Tiuj^Jl + A(uj). If the light frequency iv tends 
to infinity, then it follows from Eq.(10.20) that E — > nTiui. Since the photon momentum p = Mi^ the velocity v of the 
photon with a frequency-independent rest mass approaches c, i.e., the particle (photon) velocity v = J? — > c. Thus 
we show that in the artificial composite metamaterials such as 2TDLM media the photon velocity does not exceed 
the speed of light in a vacuum, no matter whether the phase (and group) velocity of light inside it is a superluminal 
speed or not. 

(iii) The frequency-independent rest mass of photons is useful in discussing the plus and minus signs of phase shifts 
in light scattering inside media. The phase shifts and its signs contain the information on the scattering between wave 
and potentials. According to Eq.(10.22), if A(v) is positive (negative), then the interaction energy V(w) between 
light and media is negative (positive). For this reason, it is possible for us to determine the signs (plus or minus) of 
phase shifts in light scattering process by calculating A(u) and the frequency-independent rest mass m e ff. Now with 
the development in design and fabrication of left-handed materials, the zero-index (e ~ 0, jj, ~ 0) materials receive 
attention in materials science and applied electromagnetism [44]. For this kind of media with the optical refractive 
index n ~ 0, from Eq. (10.20) it follows that 

AH ~ ^ ~ 1, (10.28) 

2 4 

i.e., A(uj) is approximately equal to "^"2 1 an d we can determine the sign by comparing the light frequency with 

2 

the photon effective rest mass m c ff. This, therefore, means that if the frequency uj is larger than m °^ fC , then A < 

and consequently the phase shift of light wave due to scattering acquires a minus sign; and if the light frequency oj is 

2 

less than m °^ fC , then A > and consequently the phase shift acquires a plus sign. 

XI. SCATTERING PROBLEM OF DOUBLE-LAYERED SPHERE CONTAINING LEFT-HANDED 

MEDIA 

The extinction properties of a sphere (single-layered) with negative permittivity and permeability is investigated by 
Ruppin. Since recently Wang and Asher developed a novel method to fabricate nanocomposite SiC>2 spheres (~ 100 
nm) containing homogeneously dispersed Ag quantum dots (2^5 nm ) [45], which may has potential applications 
to design and fabrication of photonic crystals, it is believed that the absorption and transmittance of double-layered 
sphere deserves consideration. In this paper, both the Mie coefficients and the electromagnetic field distributions in 
the two-layered sphere containing left-handed media irradiated by a planar electromagnetic wave are presented. 

A. Methods 

It is well known that in the absence of electromagnetic sources, the characteristic vectors such as E, B, D, H and 
Hertz vector in the isotropic homogeneous media agree with the same differential equation 

VV • C - V x V x C + k 2 C = 0. (11.1) 

The three independent vector solutions to the above equation is [46] 

L = VV>, M = VxaV>, N=^-VxM (11.2) 

K 

with a being a constant vector, where the scalar function ip satisfies V 2 V> + k 2 ip = 0. It is verified that the vector 
solutions M, N and L possess the following mathematical properties 

M = Lxa=)vxN, L • M = 0, VxL = 0, V • L = V 2 V> = -k 2 tp, V • M = 0, V • N = 0. (11.3) 
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Set M = mexp(— iuit) and N = nexp(— iojt), the vector wave functions M and N can be expressed in terms of the 
following spherical vector wave functions [46] 

m ° m " = t£^(At)P™(cos0)£># 2 ~ z„(fcr)^|p%>0i3, 

= ! ^^^(fcr)P„ m (cos0)->0i 1 + ±.-?-{rz n (kr)}-^P™ (cos 0)°?>0i 2 

777 f) 

In what follows we treat the Mie coefficients of double-layered sphere irradiated by a plane wave. 



1. Definitions 

We consider a double-layered sphere with interior radius a\ and external radius a 2 having relative permittivity 
(permeability) ei (^i) and e 2 ( respectively, placed in a medium having the relative permittivity e an d perme- 
ability /io- Suppose that the double-layered sphere irradiated by the following plane wave with the electric amplitude 
Eq along the z-direction of Cartesian coordinate system [46] 

E . „ ,., .^.„ 2n+l 



Si = a x E cxp(ik z - iut) = E exp(-iut) ^ ^ n Z + i) ( m °i« ~ in ein) . 

n=l ^ ' 

k k 00 2?i -|- 1 

Hi = a y —E exp(ik z - iut) = —E exp(-iwt) V i" " ■ r , (m^ + ) , (11.5) 



where 



c sm (7 0(7 

= - i ^- i Jn(fcbr)^(cosfl)^i + ^ [fcorj„(fcor)] ± ^— ^ [fc n„(fc r)]' P^cos^)™ 3 ^ (11.6) 

with the primes denoting differentiation with respect to their arguments and k n being ^/eo^o^- ln the region r < ai, 
the wave function is expanded as the following series 



2 1 

E t = E exp(-iwt) J2 * n J n + +l) («>iil - <n«) 

n=l ^ ' 



n=l ^ ' 



2n + 1 

i(n + 1 

00 



Mi^ frl n(n+l) 



Note that here the propagation constant in and n el ' n which have been defined in (11.6) is replaced with k\, i.e., 

V e iMi7- I n the region a x < r < a 2 , one may expand the electromagnetic wave amplitude as 

00 Oil 

E m = E exp(-i W t) £ i n ^jiy " *C + »>i - iC^il) . 

n=l ^ ' 

H m = -J^-Eo exp(-^i) £ ^"T^TTr fc^li + ia^n^ + + ia*n$ n ) . (11.8) 

Note that here in order to obtain m^ n and n^ n , the spherical Bessel functions j n in the spherical vector function 
m^ n and n^ n is replaced by the Hankel functions hn^ , namely, the explicit expressions for m^ n and can also 
be obtained from (11.6), so long as we replace the spherical Bessel functions j n in (11.6) with the Hankel functions 
hn \ Apparently, here the propagation constant in m^ n and n^ n should be replaced with fc 2 , i.e., ^/e 2 ^ 2 ^. 
In the region r > a 2 the reflected wave is written [46] 
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E r = E cxp(—iuit) 



2n + 1 



H r = --^-E cxp(-zwi) i n 



^ rc(n + l) 

n=l 

2n + 1 



(3) 

a n m oln 



6 r m (3) + ia 1 n (3) "| 

/ , - n l n i ]\ ^°n m eln + la n n oln J i 
n— 1 ^ ' 



(3) (3) 

where the propagation vector in and is replaced by fco, i.e., ^/eo^o^- 



,2. Boundary conditions 

At the boundary r = at, the boundary condition is given as follows 

ii x E m = ii x E t , ii x H m = ii x H t , 



where the unit vectors of Cartesian coordinate system agree with ii x i 2 = 13, i 2 x 13 = ii, h x ii = h- 
It follows from ii x E m = ii x E t that 

a™j n (N 2Pl ) + a*h£\N 2Pl ) = a^j n (N lPl ), 
Nib™ [N 2Pl j n (N 2Pl )]' + TViC [N 2Pl h£\N 2Pl )\ ' = N 2 tf n [N lPl j n (N lPl )}' , 

where Pl = k a 1: Ni = f^. 

In the similar manner, it follows from ii x H m = ii x H t that 

Mi {<£ [N 2Pl j n (N 2Pl )]' + a™ [N 2Pl hW(N 2Pl )]'} = P2 {a* [N lPl j n (N lPl )}'} , 

N 2Pl [b™j n (N 2Pl ) + b™h£\N 2Pl )] = N lPj2 blj n (N lPl ). 

At the boundary r — a 2 , in the same fashion, it follows from ii x (Ei + E r ) = ii x E m that 

3n ( P2 ) + a* n hV{ P2 ) = a™j n (N 2P2 ) + a™hW(N 2P2 ), 
N 2 [ P2jn ( P 2)}' + N 2 ¥ n [ P2 hiP (pa)] ' - b™ [N 2P2jn (N 2P2 )}' + 6* [n 2P2 H^ (N 2 p 2 )] ' , 

where p 2 = k a 2 , ^2 = f^- 

In the meanwhile, it follows from ii x (H; + H r ) = ii x H m that 

M2 {< [p2h^( P2 )] ' + [p 2 i„(p 2 )]'} - Mo [iV 2 P2.7„(^2P2)]' + a* [N 2P2 h£\N 2P2 )] j , 
M2&^i 1) (p 2 ) + mUpi) = N 2fM) b™j n (N 2P2 ) + N 2P ob™h£\N 2P2 ). 



(11.9) 



(11.10) 



(11.11) 



(11.12) 



(11.13) 



(11.14) 



3. Calculation of Mie coefficients 



Mie coefficient a r n : 

According to the above boundary conditions, one can arrive at the following matrix equation 



/ -Jn(N lPl ) Jn(N 2Pl ) h£\N 2Pl ) 

-V2[N lPl j n (N lPl )}' m[N 2Pl j n (N 2Pl )]' m N 2Pl h ( n' \N 2Pl ) 

3n{N 2 p 2 ) h ( n\N 2P2 ) 

\ mo [N 2P2 j n {N 2P2 )Y \N 2P2 h { n\N 2P2 ) 



- P 2 






h<h\ P 2) 

P 2h { n\ P2 ) 



\ 



< 



( 



\ 






jn{ P 2) 
\ P-2 [ P 2jn(P2)]' J 



(11.15) 
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For convenience, the 4x4 matrix on the left handed side of (11.15) is denoted by 



A = 



( An A 12 A 13 ^ 

A 21 A 22 A 23 

A 32 A 33 A 3A 

\ A i2 A i3 A 44 J 



whose determinant is 



detA = A xl A 22 A 33 A Ai - A n A 22 A 34 A 43 - A n A 32 A 23 A 44 + A n A 42 A 23 A 34 
— A 2 iAi 2 A 33 A 4 4 + + A21A32A13A44 — A 2 iA4 2 Ai 3 A 3 4. 



So, the Mie coefficient aL is of the form 



a„ = 



dctA r 
dctA 



with 



A r 



(A n 


A12 


A 13 





A21 


A 22 


A 23 








A 32 


A 33 


jn(p2) 


V 


A42 


A 43 


P-2 [P23n(P 



whose determinant is 



detA r = AnA 22 A 33 p 2 [p 2 j n (p2)]' - AnA 22 j n (p 2 )A4 3 - AnA 32 A 23 p 2 [p 2 j n (p2)]' + AnA4 2 A 23 j n (p 2 ) 
- A 21 A 12 A 33 p 2 [p 2 j n (p 2 )}' + A 21 A 12 j n (p 2 )A4 3 + A 21 A 32 A 13 p 2 [p 2 j n (p2)}' - A 21 A 42 A 13 j n (p 2 ). 



(11.16) 



(11.17) 



(11.18) 



(11.19) 



(11.20) 



Mie coefficient b r n : 

According to the above boundary conditions, one can arrive at the following matrix equation 



I -N 2 [N lPl j n {N lPl )^ N\ [N 2 p 1Jn (N 2Pl )}' Ni [N 2Pl h { n 1] (N 2Pl ) 

-iVlM2j„(iVlPl) N 2 /J,lj n (N2pl) N2Plh { n } (N 2 pi) 



V 



[N 2 p 2 j n (N 2 p 2 )}' 

N 2 p, jn(N 2 P2) 



N 2 p 2 h£\N 2P2 ) 
N 2 p Q h { ,P(N 2 p 2 ) 



-No 






P2h ( h\p2) 
-p,2hn\p 2 ) 



\ 



For convenience, the 4x4 matrix on the left handed side of (11.21) is denoted by 



B 



B\\ B12 

B21 B22 

B 32 

\ B i2 



Bi 3 
B 23 
B 33 
B43 



\ 


B 3 4 

B44 



J 



whose determinant is 
detB 





f 


\ 


b m 







K 




^2 [p2jn(P2)}' 




V 


P23n(p2) J 



B\\B 22 B 33 B44 — B\\B 22 B 3 4B4 3 — B\\B 32 B 23 B44 + B\\B4 2 B 23 B 3 4 
B 2 \B\ 2 B 33 B44 + B 2 \B\ 2 B 3 4B4 3 + B 2 iB 32 Bi 3 B 4 4 — B 2 iB 42 Bi 3 B 3 4. 



So, the Mie coefficient b r n is of the form 



bi 



dctB r 
dctB 



(11.21) 



(11.22) 



(11.23) 



(11.24) 



23 



with 



f B u B 12 B 13 \ 

B21 B22 B23 

B 32 S33 N 2 [p2jn(p2)]' 

\ B 42 B 43 P23n{p2) ) 



whose determinant is 



(11.25) 



detB r = B n B 22 B 33 p 2 j n (p 2 ) - BuB 2 2N 2 \p 2 jn{P2)}' B i3 - BnB 32 B 23 p 2 j n (p 2 ) + B n B 42 B 23 N 2 [p 2 jn(P2)]' 

- B 21 B 12 B 33 p 2 jn(p2) + B2lB 12 N 2 [p2jn(p2)]' B 43 + B 2 lB 32 B 13 p2 jn (/O2) - B2lB 42 B 13 N 2 [p2jn{P2)]' • (H-26) 



B. Reduced to the case of single-layered sphere 



In order to see whether the above Mie coefficients in the case of double-layered sphere is correct or not, we consider 
the reduction problem of double-layered case to the single-layered one when the following reduction conditions are 
satisfied: An = -A 12 , A 2 i = -A 2 2 N 4 = N 2 , pi = p 2 , Pi = P2- 

By lengthy calculation, we obtain 



detA = (A n A 23 - A 21 A 13 ) (A 42 A 34 - A 32 A i4 ) 

= {AxxA 23 - A21A13) i j„{N 2 p 2 )p 2 \p 2 h£\p 2 )] - Po [N 2 p 2 j n {N 2 p 2 )}' h£\p 2 ) 



and 



Thus 



dctA r = (A n A 23 - A 21 A 13 ) {po [N 2 p 2 j n (N 2 p 2 )]' j n (p 2 ) - j n {N 2 p 2 )p2 [ftejnifte)]'} ■ 

r _ ^0 [N2P2jn(N 2 p2)]' jn(P2) -jn(N 2 p2)P2 [pdnXf^ 

r (11 V ' (1) ' 

jn(N 2 p2)p2 P2h\ l '{p 2 ) - Mo [N 2 p 2 j n {N 2 p 2 )\ h n ' (p 2 ) 



(11.27) 



(11.28) 



(11.29) 



which is just the Mie coefficient a r n of single-layered sphere [46]. 

In the same fashion, when the reduction conditions Bu = —B\ 2 , B 2 i = —B22 N\ = N 2 , p\ = p 2 , pi = p 2 are 
satisfied, one can arrive at 



deti? — (BuB 23 — B 2 iBi 3 ) (B 42 B 34 — B 32 B 44 ) 

= (-B11-B23 - B 21 B 13 ) ([N 2 p 2 j n {N 2 p 2 )}' p 2 h n 1) { P 2) - Nlp j n (N 2 p 2 ) \p 2 h^(p 2 ) 



'} 



and 



So, 



dct£ r = {B n B 23 - B 21 B 13 ) {N 2 2 p 0jn (N 2 p2) [p2jn(p2)}' - [N 2 p 2 j n {N 2 p 2 )]' M2j„(p 2 )} • 

_ N^p j n (N 2 p 2 ) [p 2 j n (p 2 )]' - [N 2 p 2 j n (N 2 p 2 )}' p 2 jn(P2) 



[N 2 p 2 j n (N 2P2 )]' p 2 h { n\p 2 ) - Nlp jn(N2P2) [(toW (p 2 ) 

which is just the Mie coefficient ¥ n of single-layered sphere [46]. 

This, therefore, means that the Mie coefficients of double-layered sphere presented here is right. 



(11.30) 



(11.31) 



(11.32) 
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The 4x4 matrix A is 



C. All the Mie coefficients inside the double-layered sphere 

1. Mie coefficients a^, a™ and a™ 



/ A 12 A 13 \ 

A 22 A 23 

3n{P2] A 32 A 33 A 34 

\ M2 [P2jn(P2)}' A 42 A 43 A 44 J 



(11.33) 



whose determinant is 

detA = j n (p 2 )A 12 A 23 A 44 - j n (p2)A 22 A 13 A 44 - p 2 [p 2 j n (p 2 )]' A 12 A 23 A 34 + p 2 [p 2 j n (p 2 )}' A 22 A 13 A 34 . (11.34) 

So, 



t _ detA 
ln ~ dctl" 



(11.35) 



The 4x4 matrix A m is 



(An 





A 13 







A 2 i 





A 23 










3n{p2) 


A 33 


A 34 




V 


P2 [p23n{P2)\ 


' A 43 


A 44 


J 



(11.36) 



whose determinant is 

detAn = -Anj n {p 2 )A 23 A 44 + Anp.2 [p2j n (P2)]' A 23 A 34 + A 21 j n (p 2 )A 13 A 44 - A 2X p 2 [p2j n {P2)}' A 13 A 34 . (11.37) 

So, 

detAn 



a„ = 



detA 



(11.38) 



The 4x4 matrix An is 



Am — 



/An A 12 \ 

A 21 A 22 

A 32 jn(p2) A 34 

V A 42 p 2 [p 2 jn(P2)}' A 44 J 



(11.39) 



whose determinant is 

detAn = A u A 22 j n (p 2 )A 44 - A n A 22 A 34 p 2 [p 2 j n (p 2 )]' - A 21 A 12 j n (p 2 )A 44 + A 21 A 12 A 34 p 2 [p 2 j n (p 2 )}' . (11.40) 

So, 

m detAm 



detA 



(11.41) 



It is readily verified that under the reduction conditions An = — A2, A 2 \ = —A 22l N\ = N 2 , pi = p 2 , pi = p 2l 
one can arrive at 



detA = dctAn, a* = <C, detAn = 0, a* = 0, 
which means that the above Mie coefficient can be reduced to the those of single-layered sphere. 



(11.42) 
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2. Mie coefficients b x n , 6™ and 6™ 



The 4x4 matrix B t is 



B t 



( Bia Bia \ 

i?22 B23 

^2 [p2jn(P2)]' -B 32 B33 -B34 

V Vl3n{Pl) B 42 B43 -B44 / 



(11.43) 



whose determinant is 

detS t = N 2 [P2jn{p2)]' B 12 B 23 B 44 - N 2 [p2jn(P2)}' B 22 B 13 B 44 - p 2 jn(P2)B 12 B 23 B 34 + p 2 jn(p 2 )B 22 B 13 B 34 . (11.44) 

So, 

, t det.Bt 



dctB 



(11.45) 



The 4x4 matrix B m is 



B n 








B\ 3 





\ 


B21 





B23 










N 2 [P23n{p2) 


} B 33 


-B34 




V 


P-23n(p2) 


B43 


B44 


/ 



(11.46) 



whose determinant is 

detB m = -B U N 2 [p2jn(p2)}' B 23 B 44 + Bnp 2 j n (p 2 )B 23 B 34 + B 21 N 2 [p 2 jn(p2)]' B 13 B 44 - B 21 p 2 j n (p 2 )B 13 B 34 . 

(11.47) 

So, 

detB m 



dctB 



(11.48) 



The 4x4 matrix B r „ is 



/B n B 12 \ 

B21 -B22 

B32 N 2 [ P2 j n ( P2 )]' B 34 

\ B 42 p 2 jn{p 2 ) B44 j 



(11.49) 



whose determinant is 

dctB ffl = BnB 22 N 2 [p 2 j n {p 2 )}' B 44 - B xl B 22 B 34 p 2 j n {p 2 ) - B 21 B 12 N 2 [p 2 j n (P2)]' B 44 + B 21 B 12 B 34 p 2 j n (p 2 ). (11.50) 
So, 

detSm 



detB 



(11.51) 



It is readily verified that under the reduction conditions Bn = —B\ 2 , B 2 \ = —B 22 , N± — N 2l p 4 = p 2 , pi = p 2 , 
one can arrive at 



detB t = detB m , = detB^ = 0, 6* = 0, 
which means that the above Mie coefficient can be reduced to the those of single-layered sphere. 



(11.52) 
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D. Discussion of some typical cases with left-handed media involved 



In this subsection we briefly discuss several cases with left-handed media involved. 



(i) If the following conditions h n (N 2 pi) = 0, 
m 1' 

hn'(p 2 ) , N 2 = [l 2 = £2 



h ( n\N 2Pl )\ = 0, hi 1) (N 2 p 2 ) = -h£\p 2 ), [h£\N 2 p 2 ) 
— 1 and p = +1 are satisfied, then it is easily verified that 
det^4 = 0, deti? = and aL = oo, bl, = oo. 



(11.53) 



Thus in this case the extinction cross section of the two-layered sphere containing left-handed media is rather large. 



0. 



(ii) If the following conditions hn (N 2 p{) 
N 2 = p. 2 = e 2 = — 1 and ^ — +1 arc satisfied, then it is easily verified that 

detA r = 0, detB,. = and al = 0, bi = 0. 



h [ n\N 2Pl ) =0, h£\N 2P2 ) = Jn (p 2 ), h { ,P(N 2 p 2 ) = \j n (p 2 )}', 



(11.54) 



Thus in this case the extinction cross section of the two-layered sphere containing left-handed media is negligibly 
small. 

(hi) The case with N\ = p,\ = a = —1, N 2 = p 2 = e 2 = 0, N = p a = e = +1 is of physical interest, which 
deserves consideration by using the Mie coefficients presented above. 

Based on the calculation of Mie coefficients, one can treat the scattering problem of double-layered sphere irradiated 
by a plane wave. The scattering and absorption properties of double-layered sphere containing left-handed media can 
thus be discussed in detail, which are now under consideration and will be submitted elsewhere for the publication. 



XII. THREE KINDS OF COMPACT THIN SUBWAVELENGTH CAVITY RESONATORS MADE OF 
LEFT-HANDED MEDIA: RECTANGULAR, CYLINDRICAL, SPHERICAL 

Engheta suggested that a slab of metamaterial with negative electric permittivity and magnetic permeability (and 
hence negative optical refractive index) can act as a phase compensator/conjugator and, therefore, by combining such 
a slab with another slab fabricated from a conventional (ordinary) dielectric material one can, in principle, have a 1-D 
cavity resonator whose dispersion relation may not depend on the sum of thicknesses of the interior materials filling 
this cavity, but instead it depends on the ratio of these thicknesses. Namely, one can, in principle, conceptualize a 
1-D compact, subwavelength, thin cavity resonator with the total thickness far less than the conventional -| [47]. 

Engheta's idea for the 1-D compact, subwavelength, thin cavity resonator is the two-layer rectangular structure 
(the left layer of which is assumed to be a conventional lossless dielectric material with permittivity and permeability 
being positive numbers, and the right layer is taken to be a lossless metamaterial with negative permittivity and 
permeability) sandwiched between the two reflectors (e.g., two perfectly conducting plates) [47]. For the pattern of 
the 1-D subwavelength cavity resonator readers may be referred to the figures of reference [47] . Engheta showed that 
with the appropriate choice of the ratio of the thicknesses d\ to d 2 , the phase acquired by the incident wave at the left 
(entrance) interface to be the same as the phase at the right (exit) interface, essentially with no constraint on the total 
thickness of the structure. The mechanism of this effect may be understood as follows: as the planar electromagnetic 
wave exits the first slab, it enters the rectangular slab of metamaterial and finally it leaves this second slab. In the 
first slab, the direction of the Poynting vector is parallel to that of phase velocity, and in the second slab, however, 
these two vectors are antiparallel with each other. Thus the wave vector k 2 is therefore in the opposite direction of the 
wave vector ki . So the total phase difference between the front and back faces of this two-layer rectangular structure 
is k\d\ — \k 2 \d 2 [47]. Therefore, whatever phase difference is developed by traversing the first rectangular slab, it can 
be decreased and even cancelled by traversing the second slab. If the ratio of d\ and d 2 is chosen to be = then 
the total phase difference between the front and back faces of this two-layer structure becomes zero (i.e., the total 
phase difference is not 2mr, but instead of zero) [47]. As far as the properties and phenomena in the subwavelength 
cavity resonators is concerned, Tretyakov et al. investigated the evanescent modes stored in cavity resonators with 
backward- wave slabs [48] . 



A. A rectangular slab 1-D thin subwavelength cavity resonator 

To consider the 1-D wave propagation in a compact, subwavelength, thin cavity resonator, we first take into account 
a slab cavity of three-layer structure, where the regions 1 and 2 are located on the left- and right- handed sides, and 
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the plasmon-type medium (or a superconductor material) is between the regions 1 and 2. The above three-layer 
structure is assumed to be sandwiched between the two reflectors (or two perfectly conducting plates) [47]. Assume 
that the wave vector of the electromagnetic wave is parallel to the third component of Cartesian coordinate. The 
electric and magnetic fields in the region 1 (with the permittivity being ei and the permeability being /Ui) are written 
in the form 



E x i = E i sin (mkoz) , H yl 



niko 

iljJjJLl 



E i cos (mkoz) (0 < z < d\), 



(12.1) 



where fc stands for the wave vector of the electromagnetic wave under consideration in the free space, i.e., fc = ^, 
and in the region 2, where d\ + a < z < di + d 2 + a, the electric and magnetic fields are of the form 



E x2 = E 2 sin [n 2 k (z - d\ - d 2 - a)] , H, 



y2 



n 2 k 
iufj, 2 



E 02 cos [n 2 k (z — di — d 2 — a)] 



(12.2) 



where d\, d 2 and o denote the thicknesses of the regions 1, 2 and the plasmon (or superconducting) region, respectively. 
The subscripts 1 and 2 in the present paper denote the physical quantities in the regions 1 and 2. Note that here the 
optical refractive indices ri\ and n 2 arc defined to be n\ = and n 2 = ^/e 2 [i 2 . Although in the present paper we 

will consider the wave propagation in the negative refractive index media, the choice of the signs for n\ and n 2 will 
be irrelevant in the final results. So, we choose the plus signs for n\ and n 2 no matter whether the materials 1 and 
2 are of left-handcdncss or not. The choice of the solutions presented in (12.1) and (12.2) guarantees the satisfaction 
of the boundary conditions at the perfectly conducting plates at z = and z = d\ + d 2 + a. 

The electric and magnetic fields in the plasmon (or superconducting) region (with the resonant frequency being u> p ) 
take the form 



E xs = Aexp(/3z) + Bcxp(-/3z), H y 



[Aexp(Pz) - Bexp(-/3z)] 



(12.3) 



where the subscript s represents the quantities in the plasmon (or superconducting) region, and /3 
To satisfy the boundary conditions 

E x l\ z =di = £'a;s|z=(ii, E[yl\z=d! — E[ys\z=d! 

at the interface (z = d\) between the region 1 and the plasmon region, we should have 



(12.4) 



E m sin (nik di) = Aexp(/3di) + f?exp(— (3di), 



E 01 cos (nik di) — Acxp(/3rfi) — _Bexp(— (3di). (12.5) 



It follows that the parameters A and B in Eq.(12.3) are given as follows 



A = -cxp(-f3d 1 )E 01 



sin (nik di) + cos (nik di) 



B 



-exp (j3di) E i 



nifc 



0m 



cos (mfcorfi) — sin (mfcodi) 



In the similar fashion, to satisfy the boundary conditions 

E x2 \ Z —di+a E xs \ z =di+ai Ely 2 \ z —di+a Hys \ z — di+a 

at the interface (z = d\ + a) between the region 2 and the plasmon region, one should arrive at 

E02 sin (-n 2 k d 2 ) = A exp [(3(di + a)] + Bexp [-(3(di + a)] , 
n 2 k 



-E 02 cos (~n 2 k d 2 ) = A exp \j3(di + a)} — Bcxp [— (3{d\ + a)] . 



(12.6) 



(12.7) 



(12.6 



It follows that the parameters A and B in Eq.(12.3) are given as follows 



A = -^E 02 exp [-/3(rfi + a)] 



B = --exp[-(3{d 1 +a)]E 02 



sin (n 2 k d 2 ) - cos (n 2 k Q d 2 ) 



n 2 k 



P»2 



cos (n 2 k d 2 ) + sin (n 2 k d 2 ) 



(12.9) 
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Thus, according to Eq.(12.6) and (12.9), we can obtain the following conditions 



E, 



01 



sin (mkodi) + cos {n\kod\) 



01 



nikp 



0Mi 



cos (nifcodi) — sin {mkodi 



+ E 02 exp(-/3a) 



= 0, 



- i? 2 exp(/3a) 



sin (n 2 k d 2 ) - cos (n 2 k d 2 ) 

/?M2 

cos (n 2 kod 2 ) + sin (n 2 kod 2 ) = 0. 



rc 2 fcp 

/?M2 



(12.10) 



In order to have a nontrivial solution, i.e., to have £?oi and -E02 7^ 0, the determinant in Eq.(12.10) must vanish. 
Thus we obtain the following restriction condition 



[exp(/3a) + exp(-/3a)] 
+— [exp(/3a) - exp(-/3a) 



— tan (n 2 kod 2 ) H — - tan (mfeodi) 
Mi M2 



tan (nifc ai) tan (n 2 kod 2 ) + — , 







(12.11) 



for the electromagnetic wave in the three-layer-structure rectangular cavity. 

If the thickness, a, of the plasmon region is vanishing (i.e., there exists no plasmon region), then the restriction 
equation (12.11) is simplified to 



— tan (n 2 kod 2 ) + — tan (mkodi) = 0. 
Mi M2 



(12.12) 



In what follows we will demonstrate why the introduction of left-handed media will give rise to the novel design of 
the compact thin cavity resonator. If the material in region 1 is a regular medium while the one in region 2 is the 
left-handed medium, it follows that 



tan(nifcodi) — 711M2 
tan(n 2 /corf 2 ) n 2 Mi 



(12.13) 



According to Engheta [47] , this relation does not show any constraint on the sum of thicknesses of d\ and d 2 • It rather 
deals with the ratio of tangent of these thicknesses (with multiplicative constants). If we assume that u, d\ and d 2 
are chosen such that the small-argument approximation can be used for the tangent function, the above relation can 
be simplified as 



di 
d 2 



M2 

Mi' 



(12.14) 



This relation shows even more clearly how d\ and rf 2 should be related in order to have a nontrivial 1-D solution with 
frequency u for this cavity. So conceptually, what is constrained here is not d\ + d 2 . Therefore, in principle, one 
can have a thin subwavelength cavity resonator for a given frequency [47] . 

So, one of the most exciting ideas is the possibility to design the so-called compact thin subwavelength cavity 
resonators. It was shown that a pair of plane waves travelling in the system of two planar slabs positioned between 
two metal planes can satisfy the boundary conditions on the walls and on the interface between two slabs even for 
arbitrarily thin layers, provided that one of the slabs has negative material parameters [48]. 

In the following let us take account of two interesting cases: 

(i) If (3a — > 0, then the restriction equation (12.11) is simplified to 



— tan (n 2 fco<i 2 ) H — - tan (nifcodi) 



Mi 



which yields 



M2 



tan (riikodi) 



2(3 2 a 
k 



t \ 1 \ n l n 2kn 

tan (mkodi) tan (n2kod 2 ) + —z 

P 2 MiM2. 



0. 



fco ni 
P Mi' 



tan (n 2 /c d 2 ) 



fco ri2 
' P M2' 



(12.15) 



(12.16) 



If both nik di and n 2 fc d 2 are very small, then one can arrive at d\ = d 2 = from Eq.(12.16), which means 

that the thicknesses d\ and d 2 depend upon the plasmon parameter (3. 

(ii) If the resonant frequency lu p is very large (and hence /?), then it follows from Eq.(12.10) and (12.11) that 
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— tan (712^0^2) = 0, — tan (mkodi) = 0, 
Pi P2 



(12.17) 



namely, regions 1 and 2 are isolated from each other., which is a result familiar to us. 

In conclusion, as was shown by Engheta, it is possible that when one of the slab has a negative permeability, 
electromagnetic wave in two adjacent slabs bounded by two metal walls can satisfy the boundary conditions even if 
the distance between the two walls is much smaller than the wavelength [47] . 



B. A cylindrical thin subwavelength cavity resonator 

Here we will consider the restriction equation for a cylindrical cavity to be a thin subwavelength cavity resonator 
containing left-handed media. It is well known that the Hclmholtz equation V 2 E + k 2 ~E = in an axially symmetric 
cylindrical cavity (with the 2-D polar coordinates p and ip) can be rewritten as 



1 2 dE n 



V Z E V - + + k 2 E v = 0, 

p z p z dip 

V 2 E z + k 2 E z =0. (12.18) 

One can obtain the electromagnetic field distribution, E pi E v and H p , H v , in the above axially symmetric cylindrical 
cavity via Eq.(12.18). But here we will adopt another alternative way to get the solutions of electromagnetic fields in 
the cylindrical cavity. If the electromagnetic fields are time-harmonic, i.e., E(p, ip, z, t) = £(p, ip) exp[i(hz — kct)] and 
H(p, ip, z, t) = H(p, tp) exp[i(hz — kct)], then it follows from Maxwell equations that 

-ikcSp = - ( -^-^ - ihTtS) , 
e\p dip J 

~ikc£ v = - (ihHp 



e\ dp 
1 (dH v | \ n _1&H 1 
e \ dp p v p dip 



M= ^ + ^-1^ , (12.19) 



ikcHo = — \—^7r^- - ih£ u 



and 

p\p dip 

ikcHu> = — \ ih£p - ^ 
M \ dp 

ikcH z = 1 -(^ + 1 -E v - 1 -^). (12.20) 
p V d P P P op J 

Thus it is demonstrated that the electromagnetic fields £ p , £ v and H p , can be expressed in terms of £ z and H z , 
i.e., 



F -- ' -/ h^£l^^^k\ F -_! d. ld£ * fc2 ^M , 1991 n 

t p- k 2 -h 2 \ h dp + e P dp J ' tv ~ k 2 -h 2 \ n pdp e dp h [ 1 



and 

i (an, k 2 d£ z \ i ( idn z k 2 d£ z \ 

Hp -W^V\ h -di--p7 p -^)' n «-l^h 2 \ h -p^ + T^P~)- { } 

As an illustrative example, in what follows, we will consider only the TM wave (i.e., TL Z — 0) in the axi- 
ally symmetric double-layer cylindrical thin subwavelength cavity resonator. Assume that the permittivity, per- 
meability and radius of media in the interior and exterior layers of this double-layer cavity resonator are ei, 
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Ri and e 2 , p 2 , R 2 , respectively. According to the Helmholtz equation (with the boundary material be- 
ing the perfectly conducting medium, £2z\r 2 = 0), we can obtain E\ z = J m \\Jk\ — h\p \ < coamL P I an( j 

\ /I sin TYL^p I 



e 2z = [A Jm (v^r^fp) +B Nm (v^f-ifp)] I c ™2l 

the electromagnetic fields in both interior and exterior layers are of the form 



. Thus it follows from Eq.(12.21) and (12.22) that 



imh\ 



~ J n 



{kl-h\) P 
Jrn ( \Jk\- h\p 



hip 



cos mip 
sin my 



cos rrup 
sin my 

smrmp 
- cos my 



(12.23) 



H 



imk\ 



1 



{kf-hDpfj,! 



J77 



n 



ik\ 



^k\ - h\ P 



n lz = o, 



k\ - /l 2 p 



— sin rrup 
cos my 



fc 2_^2 \ { cos my 1 
1 1 / I smTO< /° J ' 



(12.24) 



and 



E 2p — 



E 2lp 

E2z 



ih 2 



imh 



1 1 2 

"■2 



AJ' m [yjkl-hlp^ + BN' m (Vfcf-^)] { S 



{kl-hl)p 



AJ„ 



k 2 



h\p + BN„ 



k\ - h\p 



\)k\-h\p ) r BX n 



k 2 



h\p 



cos rrup 
sin rrup 



cos rrup 
sin rrup 



J sin rrup 1 
[ — cosmos J ' 



(12.25) 



— 

7~t 2v = 
H2 



imkn 



1 



Of - ^2) P M2 
ik\ 1 

- ^2 ^ 2 



fc 2 - h\p 



h\p 



kl-hlp)+BN' m [Jkl-hlp 



— sin rrvp 

COS my 

cos rrup 
sin my 



0. 



(12.26) 



where J m and iV m denote Bessel functions and Norman functions, respectively, and k\ — y/eTpi^, k 2 = ^e 2 P2 7- 

By using the boundary conditions £ 2z \ p =r 2 = 0, -^l^fe = (due to the perfectly conducting medium at the 
boundary p — R 2 ), one can determine the relationship between A and B, i.e., 



B 



Jm ( \Jk\ - h\R 2 



k 2 - h\R 2 



(12.27) 



By using the boundary conditions £\ z \r\ = £ 2z \ri, £i v \ri = £ 2 <p\ri, one can obtain 

Jn 



k{ - hfRi = A, J, 



hi 



kf-tij 



J m y k\ - h\R x = 



k\ - h\R x 
h 2 



BN n 



k 2 - h\R x 



AJ m Uk 2 -h 2 R^\ + BN m Uk\-h\R x 



(12.28) 



The second equation in Eq. (12.28) is employed to determine the relation between hi and h 2 . With the help of the 
boundary conditions Hi z \ri — 'H 2 z\r\ , Hi v \ri — H2ip\ri, one can arrive at 
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y/kl - h\ Mi 



J' m Jk\ - h\R, 



n 2 



h\ M2 



AJ' m [^Jkl - hlR^j +BN' m (y/** - hlR^j 



(12.29) 



Combination of the first equation in Eq. (12.28) and (12.29), we have 



ni 



Mi V k 'f ~ h l J m (V k l - h'fRi) M2 V /fc 2 2 -^ AJ m (V^2 " t%Ri) + BN m (yk$ - /l^l) ' 



(12.30) 



which may be viewed as the restriction condition for the cylindrical cavity resonator. If we choose a typical case with 
hi — h 2 — 0, then the obtained restriction condition is simplified to be 



1 J' m (fcifli) _ 1 AJ' m (k 2 Rt) + BN' m (k 2 Ri) 



Ml Jm 

(fcii?i) p 2 AJ m {k 2 Ri) + BN m {k 2 Ri) ' 
Substitution of the relation (12.27) into (12.31) yields 

1 JL(hRi) = 1 J' m foRi) N m (k 2 R 2 ) - J m (k 2 R 2 ) N' m (kM 

Ml Jm 

(hRx) p 2 Jm (k 2 Ri) N m (k 2 R 2 ) - J m (k 2 R 2 ) N m (k 2 Ri) ' 



(12.31) 



(12.32) 



Eq.(12.32) is just the simplified restriction condition for the cylindrical cavity resonator. 

Similar to the analysis presented in Sec. I, it is readily shown that by introducing the left-handed media such 
cylindrical cavity will also act as a compact thin subwavelength resonator. The discussion on this subject will not be 
performed further in the present paper. 



C. A spherical thin subwavelength cavity resonator 



Here we will consider briefly the restriction equation for a spherical cavity to be a thin subwavelength cavity 
resonator containing left-handed media. Assume that the permittivity, permeability and radius of media in the 
interior and exterior layers of this double-layer cavity resonator are ei, Mi; Mi an d £2, M27 P2, respectively, and that 
the boundary medium at p = p 2 is the perfectly conducting material. Note that here the functions, symbols and 
quantities are adopted in the paper [49], e.g., N\, N 2 denote the refractive indices of interio and exterior layers, and 
j n and hn^ are the spherical Bessel and Hankel functions. 

According to the paper [49] , it follows from the boundary condition at p = p 2 that 



a m j n (N 2 p 2 ) + a™h£\N 2 p 2 ) = 0, 
b m [N 2 p 2 j n (N 2P2 )}' + C \N 2P2 h^{N 2P2 ) 



= 0. 



(12.33) 



The roles of Eqs. (12.33) is to determine the Mie coefficients a m and b m in terms of a™ and b™. 

At the boundary p = p\, it follows from the boundary condition ii x E m = ii x E t that one can obtain 

<Jn(N 2Pl ) + a™hW(N 2Pl ) = a*J n (N lPl ), 

NJ™ [N 2Pl j n {N 2Pl )\ + Nxb? \N 2Pl h^{N 2 p{)\ ' = N 2 b l n [N lPl j n (N lPl )]' . 



(12.34) 



The role of the first and second equations in Eq. (12.34) is to obtain the expressions for the Mie coefficients a l n and b l n 
in terms of a™ and b m , respectively. 

In the same manner, at the boundary p — pi, it follows from the boundary condition ii x H m = ii x H t that one 
can obtain 



Mi 
A^Mi 



2Pijn(N 2Pl )\ + a™ \N 2Pl h^{N 2Pl ) 



'}=M2< [N1 



Pdn(N-Lp-l)]' 



b™jn(N 2Pl ) + b™h^{N 2Pl ) = N^tfjniN^). 



(12.35) 



Insertion of the expressions for the Mie coefficients a m and b m and a* n and b l n in terms of a™ and b™ obtained by 
Eqs. (12.33) and (12.34) into Eqs. (12.35) will lead to a set of equations of the Mie coefficients a m and b m . In order 
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to have a nontrivial solutions of a™ and 6™, the determinant in Eqs. (12.35) must vanish. Thus we will obtain a 
restriction condition for the spherical cavity resonator. 

By analogy with the analysis presented in Sec. I, it is easily verified that by involving the left-handed media such 
spherical cavity will also serve as a compact thin subwavelength resonator. In a word, the possibility to satisfy the 
boundary conditions for small distances between metal plates is based on the fact that plane waves in Veselago media 
are backward waves, meaning that the phase shift due to propagation in a usual slab can be compensated by a negative 
phase shift inside a backward- wave slab [47,48]. We will not discuss further this topic in this paper. 



XIII. A QUANTUM MECHANICAL PROBLEM: CONE-ANGLE-INDEPENDENT GEOMETRIC 
PHASES AND PHOTON FREQUENCY SHIFT IN BIAXIALLY ANISOTROPIC LEFT-HANDED MEDIA 

Berry's discovery [50] that there exists a topological (geometric) phase (in addition to the dynamical phase that 
is familiar to physicists) in quantum mechanical wavefunction of an adiabatic process opens up new opportunities 
for investigating the global and topological properties of quantum evolution. Geometric phases thereby gained par- 
ticular and considerable attention of many researchers in various fields such as quantum mechanics [51], differential 
geometry [52], gravity theory [53,54], atomic and molecular physics [55-57], nuclear physics [58], quantum optics 
[59], condensed matter physics [60,61], molecular structures and chemical reaction [55] as well. More recently, many 
authors concentrated on their special attention on the potential applications of geometric phases to the geometric 
(topological) quantum computation, quantum decoherence and related topics [62,63]. It is well known that, due to its 
global and topological feature, geometric phases depends upon the solid angle subtended at the parameter space of 
Hamiltonian. For example, the adiabatic cyclic geometric phases of photons propagating inside a noncoplanar optical 
fiber is 27rer(l — cosO) with a denoting the photon helicity eigenvalue, where 2-7r(l — cos6>) is just the expression for 
the solid angle (with 29 being the cone angle) subtended by a curve traced by the direction of wave vector of light, 
at the center of photon momentum space [64,65]. In contrast, in this paper, we will propose a new geometric phase 
that is independent of cone angle, by taking into account the light propagation in a curved fiber made of biaxially 
anisotropic left-handed media. 

It should be noted that most of the recent theoretical works discussed mainly the characteristics of electromagnetic 
wave propagation through isotropic left-handed media, but up to now, the left-handed media that have been prepared 
successfully experimentally are actually anisotropic in nature, and it may be very difficult to prepare an isotropic left- 
handed medium [5,44]. Here we consider the wave propagation inside such biaxially anisotropic left-handed medium 
the permittivity and permeability tensors of which are written as follows 

/ e \ / -n \ 

(e) ik = -e , (A)ifc = p . (13.1) 

V o o e 3 y v ^3 / 

If the propagation vector of time-harmonic electromagnetic wave is k = (0, 0, k), then according to the Maxwellian 
Equations, one can arrive at k x E = (—kE 2 ,kEi,0), [({i)ikHk} = {—^Hi,fj,H 2 ,0). It follows from the Faraday's 
electromagnetic induction law V x E = — ^ that Hi = ^ and H 2 = Thus, the third component of Poynting 

vector of this time-harmonic wave is obtained 

S 3 = E 1 H 2 -E 2 H 1 = ^—(E 2 1 -El), (13.2) 

which implies that the Poynting vectors corresponding to the E\- and E 2 - fields are of the form 

S (D = jL kj S (2)^__^L k; (13.3) 

respectively. It is apparently seen from (13.3) that the direction of wave vector S^ 2 ) is opposite to that of S^. This, 
therefore, means that if \i > 0, then for the E\ field, this biaxially anisotropic medium characterized by (13.1) is like 
a right-handed material (regular material) whereas for the E 2 field, it serves as a left-handed one. 

In view of above discussions, it is concluded that inside the above biaxially anisotropic medium, the wave vectors 
of Ei- and E 2 - fields of propagating planar wave are apposite to each other. Consider a hypothetical optical fiber 
that is made of this biaxially anisotropic left-handed medium, inside which the wave vector of E\ field propagating is 
assumed to be k(t) = fc(sin 9 cos ip, sin 9 sin ip, cos 9) . If both 9 and ip are nonvanishing, then this fiber is noncoplanarly 
curved and in consequence the geometric phases of light will arise. It is readily verified from (13.3) that the wave 
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vector of E 2 - field is — k(t) = fc(sin 9' cos ip', sin9' sin <p', cos 9') with 9' = n — 9 and p' = tp + n. Note that for the 
latter case (i.e., E 2 - field), the expression for the time-dependent coefficient in photon geometric phases changes from 



/ * <p(t) 1 — cos 9(t) dt to J Q ip(t) 1 + cos 9(t') dt (because of 9 — > it — 9 and <p ^ p + ir ). In the next subsection 

these results will be useful in calculating the cone angle independent geometric phases of circularly polarized light in 
biaxially anisotropic left-handed media. 

Note that such biaxially anisotropic medium can be fabricated by current technology [5,44]. So, in what follows 
we will consider the geometric phases inside a noncoplanarly curved optical fiber made of the biaxially anisotropic 
left-handed media mentioned above. The time-dependent Schrodinger equation that governs the time development of 
the photon wavefunction in a curved fiber (composed of regular isotropic media) is written in the form 



.d\a,k(t),n a ) _ k(i) x k(i) 



dt 



k 2 



S |cr, \s.(t),n a ) 



(13.4) 



with S being the spin operator of photon fields. By making use of the Lewis-Riesenfcld invariant theory and the 
invariant-related unitary transformation formulation [66-68], we obtain the exact particular solutions \a, k(f), n a ) = 

exp \<pf'(t) V(t) \n a ) to Eq.(13.4), where \n a ) = \a, k(f = 0), n CT ) is the initial polarized photon state, V(t) = 

exp[/3(i)S+-/3*(i)S_] with (3(t) = -[6(t)/2] exp[-iip(t)], /3*(t) = -[0(t)/2] exp[itp(t)] [69,32] and S± = Si±iS 2 . The 
noncyclic nonadiabatic geometric phase is given as follows 



^\t) = y* p{t)[l - cos0(t')]dt'} (n a \ S 3 \n a ) . 



(13.5) 



Note here that the photon states \n a ) corresponding to the right- and left- handed polarized light with the photon 



occupation numbers being ur and are respectively defined to be \ur) 



\0 R ) and \n L ) 



(4) nj 



|0 L ). The 



creation operators of left- and right- handed circularly polarized light are constructed in terms of a\ and a\, i.e., 

respectively [70] , where a\ and a\ are the creation operators corresponding to the two 



and a 
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mutually perpendicular polarization vectors, which are also orthogonal to the wave vector k of the time harmonic 
electromagnetic wave. 

Note, however, that according to Eq.(13.3), in the above biaxially anisotropic left-handed media the wave vector of 
light corresponding to the two mutually perpendicular polarization vectors are antiparallel to each other. Since the 
wave vector of £a-field in such anisotropic left-handed media is antiparallel to that of .E^-field, we should first calculate 
the following expectation value (nR\a\ai\nR) , (ur \ a 2 02! tir), (wi,|aJoi|ni,) and (nz,|ci 2 a 2 |ni,) in order to obtain the 
expressions for geometric phases of left- and right- handed circularly polarized light in this peculiar biaxially anisotropic 



(n R - 1|, a 2 \n R ) 



\n R — 1) and 



left-handed medium. With the help of a\\riR) = -\/lFl n ft — 1)' ( n fll a i = 
(nR\a\ = —i\p^f(nR — 1|, we obtain (n R \a\ai\n R ) = ^ and (nR\a\a 2 \nR) = Hence the nonadiabatic noncyclic 
geometric phases of right-handed polarized photons corresponding to E\- and E 2 - fields are 



ip(t)(t ) 1 - cos 9(t 



dt 



4 2) (*) = ^ 



respectively, and their sum is 



(p{t)(t ) 1 + cos 9(t 



dt 



(13.6) 



Jfl(i) = {t) +4> { R\t) = UR f <p(t)(t)dt, 

Jo 



which is independent of the cone angle 9(t) of photon momentum k space. 
In the same fashion, we obtain a\\ni) = \PW\ n L — 1), (tt-lWi = 



(13.7) 



1) 



[n L \a 2 



^(ul - 1|, (n L \a\ai\n L ) = and a 2 \n L ) = 
^^f(nL — 1|, (nL\a\a 2 \nL) = ^j- ■ Hence the nonadiabatic noncyclic geometric phases of 



left-handed polarized photons corresponding to E\- and E 2 - fields are 
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p(t)(t ) 1 - cos 9(t 



dt 



<p{t){t ) 1 + cos 9(t ) 



d* 



respectively, and their sum is 
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(13.8) 



(13.9) 
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which is also independent of the cone angle 0(t). 

Thus the total geometric phases of left- and right- handed polarized photons is given by 

4ot(t)=Mt) + Mt) = (n R -n L ) [ <p{t)(t')dt, (13.10) 

Jo 

which differs from the total geometric phases of circularly light in the regular curved fiber only by a cone angle 9{t) 
of photon momentum k space. 

In the Chiao-Wu-Tomita fiber experiment [64,65], the geometric phases of photons is related close to the the 
geometric nature of the pathway (expressed by 9 and <p(t)) along which quantum systems evolve. However, in our 
case where the light propagating inside a biaxially anisotropic left-handed media, the geometric phase (13.10) is 
related only to the precessional frequency ip of wave propagation on the helicoid inside the curved fiber. 

Consider the Chiao-Wu adiabatic cyclic case in which the rotating angular frequency (i.e., the precessional fre- 
quency) of photon moving on the helicoid reads ip — ft with ft = -j== = [32] where d and a respectively denote 

the pitch length and the radius of the helix, and c is the speed of light. It follows from (13.10) that the photon 
frequency shift, A, is crfi, namely, for the right-handed circularly polarized photon, the frequency shift A is ft (due to 
a = +1), while for the left-handed circularly polarized photon, A is — £1 (due to a = —1). Generally speaking, one can 
obtain the photon frequency shift via some nonlinear optical effects such as Raman effect and left-right coupling of 
circularly polarized light in biaxially gyrotropic left-handed media [71], which may all be viewed as dynamical scheme. 
However, here we propose a so-called geometric scheme to achieve photon frequency shift, since the precessional fre- 
quency ft depends on both the pitch length d and the radius a of the fiber helix. This, therefore, means that this 
photon frequency shift is controllable by manipulating the spatial shape and helix radius of the curved optical fiber. 

For the present, physicists' control over the behavior of photons has spread to include the photon number, phase [72] 
and polarization [73] of light wave. If we could engineer all the degrees of freedom of photons, our technology would 
benefit. Already, the fiber communication, which simply guides light, has revolutionized the telecommunications 
industry. Apparently, it is of essential significance to control and utilize the degrees of freedom of photons (photon 
number, polarization, helicity, geometric phase, etc.) in information science and technology. In this paper, we suggest 
a geometric scheme, which can control the photon frequency shift by altering the curvature (and pitch length and 
helix radius) of the noncoplarnarly curved optical fiber. 

XIV. EXTRA PHASES OF LIGHT AT THE INTERFACES BETWEEN LEFT- AND RIGHT- HANDED 

MEDIA 

In this section, we will consider the effects of light appearing at the interfaces between left- and right- handed 
media. In order to treat this problem conveniently, we study the wave propagation inside an optical fiber which is 
periodically modulated by altering regular and negative media. Although it is doubtful whether such periodically 
modulated fibers could be designed and realized or not in experiments at the optical scale, it could be argued that 
the work presented here can be considered only a speculative one. But the method and results obtained via the use 
of this optical fiber system composed by such sequences of right- and left- handed materials can also be applied to 
the light propagation at the interfaces between left- and right- handed media in arbitrary geometric shapes of optical 
materials. In this periodical optical fiber, helicity inversion (or the transitions between helicity states) of photons may 
be easily caused by the interaction of light field with media near both sides of the interfaces between LRH materials. 
Since photon helicity inversion at the interfaces mentioned above is a time- dependent process, this new geometric 
phase arises during the light propagates through the interfaces (in the following we will call them the LRH interfaces) 
between left- and right- handed media, where the anomalous refraction occurs when the incident lightwave travels to 
the LRH interfaces, we think that, in the literature it gets less attention than it deserves. In what follows we calculate 
the photon wavefunction and corresponding extra phases (including the geometric phases) in this physical process, 
and emphasize that we should attach importance to this geometric phases when considering the wave propagation 
near the LRH interfaces. 

A. Model Hamiltonian 

We now treat the helicity reversal problem of light wave adjacent to the interfaces of left- and right- handed media. 
For convenience, let us consider a hypothetical optical fiber that is fabricated periodically from both left- and right- 
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handed (LRH) media with the optical refractive indices being — n and n, respectively. Thus the wave vector of photon 
moving along the fiber is respectively — in left-handed (LH) section and in right-handed (RH) section, where 
u and c respectively denote the frequency and the speed of light in a vacuum. For simplicity, we assume that the 
periodical length, b, of LH is equal to that of RH in the fiber. If the eigenvalue of photon helicity is a in right- 
handed sections, then, according to the definition of helicity, h = ^ • J with J denoting the total angular momentum 
of the photon, the eigenvalue of helicity acquires a minus sign in left-handed sections. We assume that at t = 
the light propagates in the right-handed section and the initial eigenvalue of photon helicity is a. So, in the wave 
propagation inside the LRH-periodical optical fiber, the helicity eigenvalue of h is then (— ) m a with m = [^|] , where 
\.nb\ re P resen ts the integer part of It is clearly seen that (— ) m stands for the switching on and off of the helicity 
reversal, i.e., the positive and negative value of (— ) m alternate in different time intervals. This, therefore, means that 
if 2fc(2*) < t < (2ft + 1) (2*), then (-) m = +1, and if (2k + 1) (2*) < t < (2fc + 2)(^), then (-) m = -1, where 
k is zero or a positive integer. It follows that the incidence of lightwave on the LRH interfaces in the fiber gives rise 
to the transitions between the photon helicity states (|+) and |— )). This enables us to construct a time-dependent 
effective Hamiltonian 

H(t) = ^(t)(S + + S.) (14.1) 

in terms of |+) and |— ) to describes this instantaneous transition process of helicity states at the LRH interfaces, 
where S+ = |+) (— | , S- = |— ) (+| ,^3 = \ (|+) (+| - |— ) (— |) satisfying the following SU(2) Lie algebraic commuting 
relations .?_] = 2 S3 and [S3,S±] = ±S±. The time-dependent frequency parameter w (t) may be taken to 
be u (t) = <;4ip(t), where p(t) = (— ) m with m = , and <^ is the coupling coefficient, which can, in principle, 
be determined by the physical mechanism of interaction between light fields and media. Since p (t) is a periodical 
function, by using the analytical continuation procedure, it can be rewritten as the following linear combinations of 
analytical functions 
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In what follows, we solve the time-dependent Schrodinger equation (in the unit U = 1) 

H(t) |*„(t)> (14.3) 

governing the propagation of light in the LRH- periodical fiber. According to the Lewis-Riesenfeld invariant theory 
[74], the exact particular solution (t)) of the time-dependent Schrodinger equation (14.3) is different from the 
eigenstate of the invariant I(t) only by a time-dependent c- number factor exp [j<f>o- (t)] , where 



4* (t) 



j\^(t')\[H(t')-i-^}\^(t'))dt' (14.4) 



with \& a (t)) being the eigenstate of the invariant J(i) (corresponding to the particular eigenvalue a) and satisfying 
the eigenvalue equation I (t) |$ CT (t)) — a |$ CT (t)), where the eigenvalue <r of the invariant I(t) is time-independent. 
Thus we have 



(*)> = exp 



\<t>* (t) 



\<t>„(t)). (14.5) 



In order to obtain |^ CT (t)}, we should first obtain the eigenstate |<j> CT (t)) of the invariant I(t). 



B. Photon geometric phases due to helicity inversions inside a periodical fiber made of left-handed media 

Here we investigate the time evolution of photon wavefunctions and extra phases due to photon helicity inversion at 
the LRH interfaces. As has been stated above, for convenience, we consider the wave propagation inside a hypothetical 
optical fiber which is composed periodically of left- and right- handed media. Now we solve the time-dependent 
Schrodinger equation. In accordance with the Lewis-Riesenfeld theory, the invariant I(t) is a conserved operator (i.e., 
it possesses time-independent eigenvalues) and agrees with the following Liouville-Von Neumann equation 
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^l + l[I(t),H(t)}=0. (14.6) 
It follows from Eq.(14.6) that the invariant I(t) may also be constructed in terms of S± and 53, i.e., 

I(t) = 2 I^sm9 (t) exp [-iif(t)]S+ + ^sin6(t)exp[i<p(t)]S- + cos<9(t)S 3 j ■ (14.7) 

Inserting Eq.(14.1) and (14.7) into Eq.(14.6), one can arrive at a set of auxiliary equations 

exp [—iip] (o cos 9 — iip sin 9^j — iuj cos 9 = 0, 

9 + u>smip = 0, (14.8) 

which are used to determine the time-dependent parameters, 9 (t) and <p (t), of the invariant I(t) [74]. 

It should be noted that we cannot easily solve the eigenvalue equation I (t) \<$>o- (t)) = a\$ a (t)}, for the time- 
dependent parameters 9 (t) and ip (t) are involved in the invariant (14.7). If, however, we could find (or construct) 
a unitary transformation operator V(t) to make (t)I{t)V{t) be time-independent, then the eigenvalue equation 
problem of I(t) is therefore easily resolved. According to our experience for utilizing the invariant-related unitary 
transformation formulation [54] , we suggest a following unitary transformation operator 

V(t)=exp\j3(t)S + -0*{t)S-], (14.9) 

where 0{t) and /3*(t) will be determined by calculating Iy = (t)I(t)V(t) in what follows. 
Calculation of Iy = V* (t) I (t)V (t) yields 

Iy =V*(t)I(t)V{t) = 2S 3 , (14.10) 

if and 0* are chosen to be (t) = — ^p- exp [—i<p (i)], 0* (t) = — exp [iip (t)]. This, therefore, means that we can 
change the time- dependent I(t) into a time-independent Iy, and the result is Iy — 2S3. Thus, the eigenvalue equation of 
Iy is Iy \a) = a \a) with a = ±1, and consequently the eigenvalue equation of I(t) is written I (t) V (t) \a) = aV (t) \a). 
So, we obtain the eigenstate \$ a (t)) of I(t), i.e., |$ CT (t)) = V (t) \a). 
Correspondingly H{t) is transformed into 



Hy (t) = V* (t) 



V(t) (14.11) 



and the time-dependent Schrodinger equation (14.3) is rewritten 

Hy(t)\* a (t))y=i?- t \*„(t))y (14.12) 

under the unitary transformation V(t), where \%l„ (t)) v = (t) (t)). 

Further analysis shows that the exact particular solution (t))y of the time-dependent Schrodinger equation 
(14.12) is different from the eigenstate \o) of the time-independent invariant Iy only by a time-dependent c- number 

factor exp [\(f>a (*)] [74], which is now rewritten as exp j J * (a\ Hy (i') \a) dt' j. 

By using the auxiliary equations (14.8), the Glauber formula and the Baker-Campbell-Hausdorff formula, it is 
verified that Hy(t) depends only on the operator S3, i.e., 

Hy (t) = {lo (t) sin 9 (t) cos ip (t) + ip (t) [1 - cos 9 (t)}} S 3 (14.13) 

and the time-dependent c- number factor exp \_j4>a (*)] is therefore exp ji (t) + (t) j, where the dynamical 
phase is 

<pW(t)=a[ uj(t') sm9 {t') cos <p{t')dt' (14.14) 
Jo 

and the geometric phase is 
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4>i s) (t) = a [ <p (t 7 ) [1 - cos 9 {t')] dt'. 
Jo 



(14.15) 



Hence the particular exact solution of the time-dependent Schrodingcr equation (14.3) corresponding to the partic- 
ular eigenvalue, a, of the invariant I(t) is of the form 



!*„(*)> =<*H J <j>™ (t) + <t>™ (t) \V{t)\a) 



(14.16) 



It follows from the obtained expression (14.15) for geometric phase of photons that, if the frequency parameter u> 
is small (i.e., the adiabatic quantum process) and then according to the auxiliary equations (14.8), 9 ~ 0, the Berry 
phase (adiabatic geometric phase ) in a cycle (i.e., one round trip, T ~ — ) of parameter space of invariant I{t) is 



(T) = 27tct(1-cos0), 



(14.17) 



where 27r(l — cos0) is a solid angle over the parameter space of the invariant I(t), which means that the geometric 
phase is related only to the geometric nature of the pathway along which quantum systems evolve. Expression (14.17) 
is analogous to the magnetic flux produced by a monopole of strength a existing at the origin of the parameter space. 
This, therefore, implies that geometric phases differ from dynamical phases and involve the global and topological 
properties of the time evolution of quantum systems. 



C. In biaxially anisotropic left-handed media 

Note that in the previous subsection, we treat the helicity reversals of single photon in electromagnetic media made 
of isotropic left- and right- handed materials. Now we consider this problem in a system composed by a sequence of 
right-handed (isotropic) and biaxially anisotropic left-handed media, the permittivity and permeability of the latter 
is given 





(e)ifc = -e , (fi)ik = n . (14.18) 



It has been verified in Sec. Ill that for the E\- field this biaxially anisotropic medium can be regarded as a right- 
handed material while for the E%- field it can be considered a left-handed one. This, therefore, implies that only the 
Ei- field propagating through this biaxially anisotropic left-handed medium will acquire an additional phase due to 
its helicity inversion at the LRH interfaces. But since the E 2 - field is not the eigenmode of the photon helicity, we 
cannot obtain the extra phases immediately by using the formula (14.4). According to the treatment in Sec. Ill, one 
can arrive at the expression for the total extra phases of circularly polarized light (with the occupation numbers of 
polarized photons being til and n R ) propagating through the system made of a sequence of isotropic regular media 
and biaxially anisotropic left-handed one, and the result is written as 

M) = [<t> {d) (t) + <f><*)(tj\ (14.19) 

with cf)<- d \t) + 0fe)(t) = f* {uj (t r ) sin0(t') cos ip (?) + <p (?) [1 - cos0 (t')]} dt'. It follows that if n R = n L , the added 
phase due to photon helicity reversal on the interfaces between left- and right- handed media vanishes. 
Based on the restriction [80] 

n 4 {2uj + Q l )(2lu + Q R )Q L Q R + ( 4 {uj + Q l ) 2 (lu + Q R ) 2 = (14.20) 

imposed on the frequency shifts of circularly polarized light, we now consider a possibility that only one of the polarized 
light, say E R , can be propagated in biaxially gyrotropic left-handed media. If the frequency shift fl L of left-handed 
polarized light is — cj — iT, i.e., uo + Ql = —iT, where T is a positive real number, then according to Eq.(14.20), one 
can arrive at 

n 4 (r 2 + u} 2 )(2us + Q R )Q R + ( 4 T 2 (lu + Q R ) 2 = 0. (14.21) 

It is easy to obtain Q R from Eq.(14.21), and the result is given Q R = —oj + 73(1 + j)~iu} with 7 = " ^aj^r - ■ Thus, 
the frequency of right-handed polarized light is 
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LJ + Q R = 72 (1 + 7) ilV, 



(14.22) 



which means that the frequency of Er is modified by a factor 72 (1 + 7) 2 . If no L-R coupling exists, i.e., £ = 0, 
then 7 tends to infinity and the factor 72 (1 +7)^2 approaches unity, and then the frequency shift fl^ of right-handed 
polarized light is vanishing, which can be easily seen in the expression (14.22). 

Note that in the case discussed above, the left-handed polarized light in this type of media exponentially decreases 
(due to the imaginary frequency u> + O^, which is —iT) while the right-handed one can be propagated, i.e., only one 
wave can be present in this media. So, in this case the additional phase acquired by photon wavefunction due to 
hclicity inversions on LRH interfaces is 

(14.23) 

Additionally, it is of interest to show that the above scheme is applicable to the detection of quantum-vacuum 
geometric phases [32]. Since the left-handed polarized light (including the quantum vacuum fluctuation correspond- 
ing to the left-handed polarized light) cannot be propagated in this media [6,16], the quantum- vacuum geometric 
phase of right-handed polarized light will not be cancelled by that of left-handed one, namely, the only retained 
quantum- vacuum geometric phase is that of right-handed circularly polarized light and therefore it is possible for the 
nonvanishing quantum-vacuum geometric phases to be detected in experiments. 

D. Discussion: physical significance and potential applications 

It is worthwhile to point out that the geometric phase of photons caused by helicity inversion presented here is of 
quantum level. However, whether the Chiao-Wu geometric phase due to the spatial geometric shape of fiber is of 
quantum level or not is not apparent (see, for example, the arguments between Haldane and Chiao et al. about this 
problem [77,79]), since the expression for the Chiao-Wu geometric phase can be derived by using both the classical 
Maxwell's electromagnetic theory, differential geometry and quantum mechanics [78,77,79]. However, the geometric 
phase in this paper can be considered only by Berry's adiabatic quantum theory and Lewis-Riescnfcld invariant theory, 
namely, the classical electrodynamics cannot predict this geometric phase. Although many investigators have taken 
into account the boundary condition problem and anomalous refraction in left-handed media by using the classical 
Maxwell's theory [7,16], less attention is paid to this geometric phase due to hclicity inversion. It is believed that 
this geometric phase originates at the quantum level, but survives the correspondence-principle limit into the classical 
level. So, We emphasize that it may be essential to take into consideration this geometric phase in investigating the 
anomalous refraction at the LRH interfaces. 

It is well known that geometric phases arise only in time-dependent quantum systems. In the present problem, the 
transitions between helicity states on the LRH interfaces, which is a time-dependent process, results in the geometric 
phase of photons. This may be viewed from two aspects: (i) it is apparently seen in Eq.(14.8) that if the frequency 
parameter u> in the Hamiltonian (14.1) vanishes, then <p = and the geometric phase (14.15) is therefore vanishing; 
(ii) it follows from (14.2) that the frequency coefficient u (t) of Hamiltonian (14.1) is 

-&£['-<->*] «(£«)■ <14 - 24 » 

k=l 

Since |cos (^T^) | < 1, the frequency coefficient, the transition rates between helicity states, and the consequent time- 
dependent phase (tp^ (i) + (t)) greatly decrease correspondingly as the periodical optical path nb increases. Thus 
we can conclude that the interaction of light fields with media near the LRH interfaces gives rise to this topological 
quantum phase. 

In addition to obtaining the expression (14.15) for geometric phase, we obtain the wavefunction (14.16) of photons 
in the LRH- optical fiber by solving the time-dependent Schrodinger equation (14.3) based on the Lewis-Riesenfcld 
invariant theory [74] and the invariant-related unitary transformation formulation [75,76]. We believe that this would 
enable us to consider the propagation of light fields inside the optical fiber in more detail. 

In the above treatment, we constructed an effective Hamiltonian (14.1) to describe the time evolution of helicity 
states of photons. It should be noted that the method presented here is only a phenomenological description of 
propagation of lightwave in the LRH-periodical fiber. This phenomenological description is based on the assumption 
that the direction of wave vector k becomes opposite nearly instantaneously on the LRH interfaces. This assumption 
holds true so long as the periodical length b is much larger than the wavelength of lightwave in the fiber. 
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XV. WAVE PROPAGATION IN GENERALIZED GYROTROPIC MEDIA 



In this section, we investigate the wave propagation in a generalized gyrotropic medium with the following permit- 
tivity and permeability 
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(15.1) 



This type of materials is a generalization of the ordinary gyrotropic medium with 
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(15.2) 



and can encompass the electromagnetic parameters of and Mife m (15-1) material responses experimentally obtained. 

Assume that the wave vector k of electromagnetic wave propagating inside this medium is parallel to the the third 
component of the Cartesian coordinate system. According to Maxwcllian Equations, one can arrive at 



d 2 d 2 

-V 2 E 1 = fl ( M21^£>2 - M22^£>1 



d 2 d 2 
-V 2 E 2 = -no [ m-gpD 2 - W-QpDi ) , 



where the electric displacement vector D\ and D 2 are of the form 

D\ = e (e\\E\ + ei 2 E 2 ) , D 2 = e (e 2 iEi + e 22 E 2 ) . 
Due to the transverse nature of planar electromagnetic wave, the following acquirements are satisfied 

E 3 = 0, £f 3 =0, k • E = 0, k • H = 0, k • D = 0, k • B = 0. 
Thus with the help of above equations, it is verified that 



\ 1 E x = -eoMo 



V 2 E 2 = e Mo 



d 2 d 2 ' 

(M21 e 21 — M22£ll) ^2^! + (M21 e 22 — M22£l2) ^2^ 2 

d 2 d 2 

(Mll e 21 — Ml2 £ ll) + (Mll e 22 — Ml2 £ 12) ~QfiE 2 



So, 



2 / Ei ± iE 2 



1 d 2 

^ j = ^eoMo [(M22£h - M2l£2l) ± « (Mll e 21 - Ml2£ll)] 

1 d 2 

+ ^| e 0M0 [(M22ei2 - M2l£22) ± % (Mll e 22 - ^Yl^\l)\ E 2 . 



It follows from (15.1) that 



where 



(M22E11 - M2l£2l) + i (mh £ 2i - Mi2£n) = A + + iB + , 

(M22C12 - M21C22) + % (Mll e 22 - Ml2£l2) = iA+ + B + , 



A + = a c ah — a e Ch - d c dh — id c bh — b c bh + ib e dh + c Ch — c c ah, 
B + = -a c 6 h - id c a h + id c c h - b c c h + b c a h - c 6 h + ic c d h + ia e d^. 



In the same fashion, one can arrive at 



(M22£ii - M2i£2i) - i (Mn e 2i - Mi2£n) = A_+ iB_ 

(M22C12 - M2l£22) - i (Mll e 22 - Ml2£l2) = ~%A- - B- 



(15.3) 



(15.4) 



(15.5) 



(15.6) 



(15.7) 



(15.8) 



(15.9) 



(15.10) 
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where 



A- = a c a h + a e Ch - d c d h + id c b h - b c b h - ib c d h + c c c h + c c a h , 
B_ = a c fe h - id e dh - id e Ch - b c c h — fe a h — c c b h — ic c d h + ia e dh- 

Hence, the wave equations of left- and right- handed circularly polarized light 12 is written 



(15.11) 



V 2 
V 



Ei - iE 2 

2 f E 1 +iE 2 
V2 



V2 

cP_ ( E 1 +iE 2 
dt 2 { y/2 



_ d 2 (E l +iE 2 
dt 2 



V2 
E 1 - iE 2 



V2 



(15.13) 



It should be pointed out that the interaction between left- and right-handed polarized light arises in Eq.(15.13), 
which is analogous to the Josephson's effect and may therefore be of physical interest. 

The coupling of left-handed polarized light to the right-handed one can be characterized by the frequency shift f2 ^ 
and il R , namely, the amplitudes of left- and right- handed circularly polarized light propagating along the z-direction 
are written as follows 



Er 



<7 



(u + n L )t- y[Al-z + $L 



E R ~ cxp { - 



(uj + fi fl ) t - ^JA^-z + cf) R 



(15.14) 



Insertion of Eq.(15.14) into Eq.(15.13) yields 



A_uj 2 E l = A- (u + n L ) E L + tS_ (w + Q R ) 2 E R , A+uu 2 E R = A+ (w + Q R f E R + iB+ (w + Q L ) 2 E L , (15.15) 



and consequently 



a+a_ (2ujn R + n 2 R ) (2M L + nl) + b + B- (w + n L f {uj + n R f = o, 



(15.16) 



which is a restricted condition regarding the frequency-shift relation between left- and right- handed circularly polar- 
ized light. 

Note that in the case of conventional gyrotropic media characteristic of such permittivity and permeability tensors 
(15.2), where d c = b c = d^ = b^ = 0, a c = e\, c c = — e 2 , ah = fii, Ch = — M2, the optical refractive index squared is of 
the form 



A± = (ei ± e 2 ) (m ± fi 2 ) , 
and the coupling coefficients B± are vanishing. 



(15.17) 
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12 For simplicity, without loss of generality, it is assumed that the two mutually perpendicular real unit polarization vectors 
e(k, 1) and e(k,2) can be taken to be as follows: ei(fc, 1) = E2(k,2) = 1, ei(k, 2) = S2(k, 1) = and es(k, 1) = £3(^,2) = 0. 

dA 

at 



Thus by the aid of the formula E = — ^ for the electric field strength, in the second-quantization framework one can arrive 
at [49] 



Er = El ^J^ 2 = i J d 3 ky/^^-p-[a i/ (fc) exp(— ik ■ x) — a R (k) exp(ik ■ x)], 



E L = El ^J^ 2 — * J d^ky^^^g [ajj(fc) exp(— ik ■ x) — a L (k) exp(ik ■ x)]. (15.12) 
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Appendix A. 

In this Appendix, we will briefly consider the permittivity of the metal thin wire array medium. For a nearly free 
electron in the structures of array of long metallic wires (ALMWs) [18], the equation of motion is ^ = ^E, and the 

electric current density j therefore satisfies the equation || = ^-E, where j = new with n being the electron number 
density in ALMWs media. If the differential operator ^ can be replaced with iu, where to is the frequency of incident 
electromagnetic wave, then we can obtain 

-E. (15.18) 



lurni 



Thus with the help of j = ^ = iwP, one can arrive at 

2 



P = -2^-E. (15.19) 



According to the definition of e, i.e., P = (e — l)eoE, it follows that 

s 2 
me 



_E with ^ 2 = £— . (15.20) 



Appendix B. 

SRR essentially behaves as a LC circuit, the equation of which may be as follows 

according to the electromagnetic induction law. By replacing the derivative operator ^ with ico, and using the electric 
current I = -^q = ituq (hence q = one can arrive at 

B*r£ 2 

' = - L R . 1- ( 15 - 22 ) 



w a - ifyto 



2 _ ;R, 

L LC 



Thus we have the magnetization (total moments per unit volume) 



M = Nlirrl = — \ (15.23) 



which is equal to j^^-B. So, we may obtain 



1 1 M-g) 

= 1 it £ , (15.24) 



which is of the form 



F^ 2 
cj 2 + iuoY — Lli'q ' 
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